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We present measurements of e^e“ production at midrapidity in Au+Au collisions at .^/sfffj = 

200 GeV. The invariant yield is studied within the PHENIX detector acceptance over a wide range of 
mass (mee < 5 GeV/c^) and pair transverse momentum (pr < 5 GeV/c), for minimum bias and for 
five centrality classes. The e''"e“ yield is compared to the expectations from known sources. In the 
low-mass region (rriee = 0.30-0.76 GeV/c^) there is an enhancement that increases with centrality 
and is distributed over the entire pair pr range measured. It is significantly smaller than previously 
reported by the PHENIX experiment and amounts to 2.3±0.4(stat) ±0.4(syst) ±0.2™°'^®' or to 1.7± 
0.3(stat)±0.3(syst)±0.2™°‘^'®* for minimum bias collisions when the open heavy flavor contribution is 
calculated with pythia or mc@nlo, respectively. The inclusive mass and pt distributions as well as 
the centrality dependence are well reproduced by model calculations where the enhancement mainly 
originates from the melting of the p meson resonance as the system approaches chiral symmetry 
restoration. In the intermediate-mass region (mee = 1.2-2.8 GeV/c^), the data hint at a significant 
contribution in addition to the yield from the semileptonic decays of heavy flavor mesons. 

PACS numbers: 25.75.Dw 


I. INTRODUCTION 

Dileptons are important diagnostic tools of the Quark 
Gluon Plasma (QGP) formed in ultra-relativistic heavy 
ion collisions [1]. They are unique observables for their 
sensitivity to the chiral symmetry restoration phase tran¬ 
sition expected to take place together with, or at similar 
conditions to, the deconfinement phase transition HE]. 
When chiral symmetry is restored, the chiral doublets, 
such as the p and the ai mesons, become degenerate in 
mass. As the oi meson is very difficult to observe experi¬ 
mentally, the p meson is the main observable in this con¬ 
text. Due to its very short lifetime (t ^ 1.3 fm/c), the p 
meson quickly decays after its formation and is therefore 
a sensitive probe of the medium where it is formed. The 
p meson is mostly produced close to the phase boundary 
and possible modifications of its spectral function in the 
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high temperature and density conditions prevailing there 
are thus imprinted in its decay products. The decay into 
dileptons, as opposed to hadrons, is of particular interest 
as they escape unaffected by the interaction region, thus 
carrying this information to the detectors. 

Dileptons are sensitive to the thermal radiation emit¬ 
ted by the system, both the partonic thermal radiation 
(quark annihilation into virtual photons, gg —>■ 7 * —>■ 
l^l~) emitted in the early stage of the collisions as well 
as the thermal radiation emitted later in the collision 
by the hadronic system. The main channel of the lat¬ 
ter is pion annihilation, mediated through vector meson 
dominance by the p meson ( 7 r“'" 7 r“ —> p —> 7 * —> 
Dileptons are produced by a variety of sources all along 
the entire history of the collision and it is necessary to 
know precisely all these sources in order to single out the 
interesting signals characteristic of the QGP related to 
chiral symmetry restoration or thermal radiation [4]. 

The CERES experiment pioneered the study of di¬ 
electrons at the Super Proton Synchrotron (SPS). A 
strong enhancement of low-mass electron pairs (wee < 
1 GeV/c^) with respect to the cocktail of expected 
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hadronic sources, was found in all nuclear systems stud¬ 
ied, in S-I-Au collisions at 200 AGeV [5], in Pb-|-Au col¬ 
lisions at 158 AGeV [SI |7| and in Pb-|-Au collisions at 
40 AGeV [5]. The enhancement was confirmed and fur¬ 
ther studied by the high statistics NA60 experiment that 
measured dimuons in In-|-In collisions at 160 AGeV i- 
dl]. In both experiments, the low-mass dilepton enhance¬ 
ment is explained by in-medium modification of the p 
meson spectral function [ISIllHj . The data rule out the 
conjectured dropping mass of the p meson as the sys¬ 
tem approaches chiral symmetry restoration [MS]- In¬ 
stead, the data are well reproduced by a scenario in which 
the p meson copiously produced by tt+tt” annihilation 
is broadened by the scattering off baryons in the dense 
hadronic medium. The low-mass dilepton excess is thus 
identified as the thermal radiation signal from the hadron 
gas phase with a modified p meson spectral function. A 
recent paper shows that in-medium modifications of vec¬ 
tor and axial vector spectral functions lead to degeneracy 
of the p and ai meson masses providing a direct link be¬ 
tween the broadening of the p meson spectral function 
and the restoration of chiral symmetry |22j . 

NA60 found also an excess at higher masses {mi+i-=l- 
3 GeV/c^). Using precise vertex information this excess 
was associated with a prompt source originating at the 
vertex, as opposed to semi-leptonic decays of D mesons 
that originate at displaced vertices. The excess can be 
explained as thermal radiation from the QGP (aHnidS] 
but other interpretations based on hadronic models, sim¬ 
ilar to those that explain the low mass excess [I31II11, or 
on hadronic rates constrained by chiral symmetry con¬ 
siderations m can also reproduce the data. 

At the Relativistic Heavy Ion Collider (RHIC), the 
PHENIX experiment reported a strong enhancement 
of low mass pairs in Au-|-Au collisions at = 

200 GeV [23]. In the 0%-10% most central collisions, 
where the excess is concentrated, the enhancement fac¬ 
tor, defined as the ratio of the measured yield over the 
cocktail yield reaches an average value of 7.6±0.5(stat)± 
1.3(syst) ± 1.5 (cocktail) in the mass range niee = 0.15- 
0.75 GeV/c^. All models that successfully reproduce the 
SPS results fail to explain the PHENIX data [23[ [24] . 

The PHENIX result [23] was characterized by a con¬ 
siderable hadron contamination of the electron sample 
and by a small signal to background (S/B) ratio. In 
an effort to improve upon this measurement, a hadron- 
blind detector (HBD) was developed and installed in the 
PHENIX experiment [2H1 - [27| . The HBD provides addi¬ 
tional electron identification, additional hadron rejection 
and improves the signal sensitivity. 

In this paper we present dielectron results obtained 
with the HBD in 2010 for Au-|-Au collisions at -^/s^ 
— 200 GeV. The paper is organized as follows. Section 
[IT] describes the PHENIX detector with special emphasis 
on the HBD. In Section [Hi] we give a detailed account of 
the various steps of the data analysis including electron 
identification, pair cuts and background subtraction that 
is the crucial step in this analysis. The raw mass spec¬ 


tra, efficiency corrections and systematic uncertainties of 
the data are also discussed in this section. Section Hvl 
describes the procedures used to calculate the expected 
dielectron yield from the known hadronic sources. The 
results, including invariant mass spectra, pt distributions 
and centrality dependence, are presented in SectionJV] In 
the same section, the results are discussed with respect 
to previously published results and compared to available 
theoretical calculations. A summary is given in Section 

m 


II. PHENIX DETECTOR 

Figure[^shows a schematic beam view of the PHENIX 
central arm detector, as used during 2010 data taking. 
A detailed description of the detector, except the HBD, 
can be found in |28j . In this section, we give only a brief 
description of the PHENIX sub-systems relevant for the 
present analysis: global detectors, central magnet, cen¬ 
tral arm detectors, including drift chambers (DC), pad 
chambers (PC), ring-imaging Cerenkov (RICH) detec¬ 
tors, time-of-flight (TOE) detectors and electromagnetic 
calorimeters (EMCAL) and the HBD. 



West Beam View East 


FIG. 1. (Color online) Beam view (at a = 0) of the PHENIX 
central arm spectrometers during 2010 data taking. 


A. Global detectors 

The measurement of the collision-vertex position, time, 
and centrality, as well as the minimum-bias (MB) trig¬ 
ger, is provided by two beam-beam counters (BBC) [29] . 
Each BBC comprises 64 quartz Cerenkov counters, lo¬ 
cated at ±144 cm along the beam axis from the center of 
PHENIX, with 27r azimuthal coverage over the pseudo¬ 
rapidity interval 3.0 < \ri\ < 3.9. The collision-vertex po¬ 
sition along the beam direction z is determined from the 
difference of the average hit time of the photomultiplier 
tubes (PMTs) between the north and the south BBC. 
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The z-vertex resolution ranges from ^0.5 cm in central 
Au+Au collisions to ~2 cm in p+p collisions. The MB 
trigger requires a coincidence between at least two hits 
in each of the BBC arrays thus capturing 92 ± 3% of the 
total inelastic cross section m- 


B. Central magnet 


The PHENIX central magnet comprises two pairs of 
concentric coils, an inner coil pair and an outer coil pair, 
that can be operated independently and create an axial 
magnetic field parallel to the beam axis [31] . The coils are 
usually operated with current flowing in the same direc¬ 
tion (the or-configuration) so that their magnetic 

fields add together. For the dilepton measurement with 
the HBD in the 2010 run, the coils were operated with 
equal currents flowing in opposite directions. In this so 
called H— configuration, the inner coil counteracts the 
action of the outer coil so that their magnetic fields can¬ 
cel each other, creating an almost field free region in the 
inner space extending from the beam axis out to a radial 
distance of ^60 cm where the inner coil is located (see 
Fig. 1 of Ref. [IZj). The field free region preserves the 
opening angle of e“'"e“ pairs and this is an essential pre¬ 
requisite for the operation of the HBD. The HBD exploits 
the fact that the opening angle of e+e“ pairs originating 
from 7 conversions or from tt^ Dalitz decays is very small. 
When only one of the two tracks is reconstructed in the 
central arms, the HBD can reject them by applying an 
opening ang le cut or a double signal cut on the HBD hits 
(see Section H Dj ). In this configuration however, the to¬ 
tal field integral is J B ■ dl — 0.43 Tm, about 40% of the 
value in the -|--|- configuration. 


C. Central arm detectors 

PHENIX measurements at midrapidity are made with 
two central arm spectrometers, as shown in Fig. Each 
central arm covers pseudorapidity |?7| < 0.35 and az¬ 
imuthal angle Acj) = 7r/2. 

Charged-particle tracks are reconstructed using hit in¬ 
formation from the DC, the first layer of PC (PCI) and 
the collision point along the z-direction [35] . The DCs 
are located outside the magnetic field in the radial dis¬ 
tance 2.02-2.46 m from the beam axis. They provide 
an accurate measurement of the particle trajectory in 
the plane perpendicular to the beam axis. The PCls 
are multiwire proportional chambers located just behind 
the DC at 2.47-2.52 m in radial distance from the beam 
axis [33j . They provide a three dimensional space point 
that is used to determine the track origin along the beam 
axis. The transverse momentum { pt ) of each particle is 
determined from the bending of its trajectory in the az¬ 
imuthal direction. The total momentum p is determined 
by combining px with the polar angle information of PCI 
and the vertex position z. The reconstructed tracks are 


projected onto the HBD (see next subsection) and onto 
the central-arm detectors that provide electron identifi¬ 
cation: RICH, EMCal, and TOF. 

The RICH is the primary central-arm detector used for 
electron identification in PHENIX [Mj, and is located in 
the radial region of 2.5-4.1 m, just behind PCI. The 
RICH uses CO2 as the gas radiator at atmospheric pres¬ 
sure, and has a Cerenkov threshold of 7 = 35. This 
corresponds to a momentum threshold of 18 MeV/c for 
electrons and 4.7 GeV/c for pions. Two spherical mirrors 
reflect the Cerenkov light and focus it onto two arrays of 
1280 PMTs each located outside the acceptance on each 
side of the RICH entrance window. The average number 
of hit PMTs per electron track is ^5, and the average 
number of photo-electrons detected is ~10. Below the 
pion threshold, the pion rejection is ^10^ in p+p or low 
multiplicity collisions. However, in high-multiplicity col¬ 
lisions, hadron tracks are misidentified as electrons when 
their trajectory is nearly parallel to that of a genuine 
electron. This effect limits the c/tt separation to ^10“^ 
in central Au-f An collisions and requires special care as 
described below. 

The EMCal measures the energy deposited by elec¬ 
trons and their shower shape |35j . It comprises eight 
sectors each covering A(j) « 7r/8 in azimuth, where 
six sectors are made from lead-scintillator (PbSc) with 
an energy resolution 4.5% © 8.3%/a/£’ [GeV] and two 
are lead-glass (PbGl) with an energy resolution 4.3% © 
7.7 %/a/E [GeV]. The radial distance from the beam axis 
is 5.10 m for PbSc and 5.50 m for PbGl (see Fig. [^. The 
matching of the measured energy to the track momentum 
is used to identify electrons. The latter are all relativis¬ 
tic in the accepted momentum range {pr > 0.2 GeV/c), 
hence the energy-to-momentum ratio is close to unity. 

To further separate electrons and hadrons we use the 
time-of-flight information from the PbSc part of the EM- 
Gal which covers 75% of the acceptance but has a valid 
time response for 64% of the acceptance. In addition, 
we use the time-of-flight information from the TOF-east 
detector (TOF-E) |36| covering an additional 16% of the 
acceptance. The former has a time resolution of ^450 ps, 
while the latter has a resolution of ^150 ps. The rest of 
the acceptance, 9%, does not have a usable TOF cover¬ 
age, because the time resolution of ^700 ps provided by 
PbGl detectors is not sufficient for an effective separation 
of electrons and hadrons. 


D. The Hadron Blind Detector 

The HBD was installed in PHENIX prior to 2010. A 
detailed description of the concept, construction and per¬ 
formance of the HBD is given in Ref. HZj. Only a brief 
account is given here with emphasis on the specific as¬ 
pects relevant to the present analysis. 

The HBD provides additional electron identification 
and additional hadron rejection to the central arm de¬ 
tectors. Its main task is to recognize and reject 7 con- 
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versions and tt*^ Dalitz decays which are the dominant 
sources of the combinatorial background. Very often, 
only one of the two tracks of an e“'"e“ pair from these 
sources is detected in the central arm, whereas the sec¬ 
ond one is lost because it falls out of the acceptance, is 
curled by the magnetic field or is not detected due to the 
inability to reconstruct low momentum tracks with < 
200 MeV/c. The HBD exploits the fact that most of these 
pairs have a very small opening angle and thus produce 
two overlapping hits in the HBD, resulting in a charge 
response with an amplitude double the one correspond¬ 
ing to a single hit. Being sensitive to electrons down to 
very low momentum (see below), the HBD can detect 
both tracks and can effectively reject them by applying 
a double hit cut on the HBD signal. On the other hand, 
decays with a large opening angle between the electron 
and positron produce two well separated single hits on 
the HBD pad plane as illustrated in Fig. The ability 
to distinguish single from double hits is one of the main 
performance parameters of the HBD. This is illustrated 
in Fig. 1^ which shows the HBD response to single and 
double electron hits in real data. Single and double hits 
are selected from reconstructed low-mass pairs with large 
(> 100 mrad) and small (< 50 mrad) opening angles, re¬ 
spectively. 



FIG. 2. (Color online) Sketch illustrating the HBD response 
to an e''"e“ pair from 7r° Dalitz decay and from a cj) meson 
decay. The circles represent the Cerenkov blobs whereas the 
hexagons are the hexagonal pads of the HBD readout plane. 

The HBD is a Cerenkov detector. It has a 50 cm long 
radiator directly coupled, in a windowless configuration, 
to a triple gas-electron-multiplier (GEM) detector 
which has a Csl photocathode evaporated on the top 
face of the upper-most GEM foil and pad readout at 
the bottom of the GEM stack (see Fig. |^. The HBD 
uses pure CF4 at atmospheric pressure that has an av¬ 
erage Cerenkov threshold of 7 = 28.8 over the detector 
bandwidth, corresponding to a momentum threshold of 
~ 15 MeV/c for electrons and ^ 4.0 GeV/c for pions. In 
this scheme, Cerenkov radiation from particles passing 
through the radiator is directly collected on the photo¬ 
cathode forming a circular blob image rather than a ring 
as in a RICH detector. The pad readout plane comprises 



FIG. 3. (Color online) HBD response to single electron hits 
and double electron hits in the 60%-92% centrality bin. The 
two distributions are normalized to give an integral yield of 
one. 


hexagonal cells with a hexagon side of 1.55 cm. One cell 
subtends an opening angle of approximately 50 mrad and 
has an area of 6.2 cm^, comparable to the blob size which 
has a maximum area of 10 cm^. The electron response of 
the HBD is thus typically distributed over a maximum of 
3 readout cells and subtends a maximum opening angle 
of 75 mrad. 



FIG. 4. (Color online) Triple GEM stack operated in reverse 
bias mode where ionization electrons produced by a charged 
particle are repelled toward the mesh. 

The hadron blindness property of the HBD is achieved 
by operating the detector in reverse bias mode where 
the mesh defining the detection volume is set at a lower 
voltage with respect to the Csl photocathode [23111] (see 
Fig. 1^. Consequently, the ionization electrons produced 
by charged particles in the drift region defined by the 
entrance mesh and the photocathode are mostly repelled 
towards the mesh. Only the ionization electrons created 
in a thin layer of ^100 fim above the photocathode are 
collected and amplified by the GEM stack leading to a 
very small signal, equivalent to a few p.e., localized in 
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one single cell of the pad plane. 

The choice of CF 4 in a windowless configuration as the 
common gas for the radiator and the detector amplifica¬ 
tion medium, results in a large bandwidth of UV photon 
sensitivity from 6.2 eV (the threshold of the Csl photo¬ 
cathode) up to 11.1 eV (the CF 4 cut-off). This trans¬ 
lates into an average yield of 20 photo-electrons (p.e.) 
per electron, as shown in Fig. corresponding to a mea¬ 
sured figure of merit Nq of 330 cm“ ^, very high for a gas 
Cerenkov detector m- 

The HBD is located close to the interaction vertex, 
in the field-free region, starting immediately after the 
beam pipe at r = 5 cm and extending up to r = 60 cm. 
The detector comprises two identical arms, each covering 
112.5° in azimuth and ±0.45 units of pseudorapidity. The 
active area of each arm is subdivided into 10 detector 
modules, 5 along the azimuthal axis and 2 along the 2 ; 
axis. With this segmentation, each detector module is ^ 
23 X 27 cm^ in size. The material budget (See Table|^ in 
front of the GEM detectors is 0.62% of a radiation length 
dominated by the CF 4 contribution of 0.56%. To this, 
one has to add the contribution of the GEM stack, the 
vessel back plane and the front-end electronics attached 
to the vessel to give a total of 2.4% of a radiation length 
for the entire detector. 


TABLE I. Material budget of the HBD within the central arm 
acceptance m- 


Component 

Radiation length 

(%) 

Window (aclar/kapton) 

0.04 

Gas (CF 4 ) 

0.56 

GEM stack 

0.42 

Vessel back plane ± front-end electronics 

1.4 

Total 

2.4 


Good gain calibration is crucial to achieve the best 
possible separation between single and double hits in the 
HBD. Gain variations occur as a function of time due to 
two main factors: (i) variations of temperature and pres¬ 
sure and (ii) charging effects of the GEM foils produce an 
initial rise of the gain after switching on the HV, that can 
last for several hours before stabilizing [35]. These gain 
variations are taken into account by performing a gain 
calibration of each module every three minutes during 
data collection. This is done by exploiting the scintilla¬ 
tion light produced by charged particles traversing the 
CF 4 radiator. The scintillation signal is easily identified 
by the characteristic exponential shape of single electrons 
in the HBD pulse height distribution of low-multiplicity 
Au±Au collisions m- Furthermore, the average cell 
charge per event was found to slowly decrease by 10 %- 
15% over the 10 week duration of the run for some of the 
modules. This is attributed to a slow deterioration of 
the quantum efficiency of the photocathodes. This effect 


was noticed in ^40% of the modules, the others did not 
show any sign of aging although all photocathodes were 
produced under identical procedures. An additional time 
dependent correction factor is applied to account for this 
effect. 

In high multiplicity Au±Au collisions, a large amount 
of scintillation light is produced by charged particles 
traversing the GF 4 gas, resulting in a large detector occu¬ 
pancy. The number of photoelectrons per cell can be as 
high as ~10 in the most central collisions. This underly¬ 
ing event background is subtracted on an event-by-event 
basis. For each event and for each module the average 
charge per unit area (Q) is calculated as: 

w) = E Qcell/ ^cell 1 (f) 

where Qceii and Gceii are the cell charge and area, respec¬ 
tively. The summation is carried out over all the cells of a 
given module, excluding the cells that are matched to an 
electron track and their first neighbors. The cell charge 
used for further analysis Q^eih then given by: 

Qcell ~ Qcell {Q) X CLcell (2) 

After subtraction of the underlying event charge, two 
independent algorithms are used for the HBD hit recog¬ 
nition. The first is a stand-alone algorithm in which a 
cluster is formed by a seed cell with Q*^ii > 3 p.e. to¬ 
gether with the fired cells (defined as > 1 p.e.) 

among its first six neighbors. Such clusters can have up 
to seven cells. A central arm electron track projected 
onto the HBD readout plane is then matched to the clos¬ 
est cluster. This algorithm works very well in p+p or 
peripheral Au±Au collisions producing a typical single 
electron response with an average of 20 p.e.. In higher 
multiplicity events, this algorithm yields a higher charge 
per electron and a higher fraction of fake hits as it picks 
up more charge from the fluctuations of the underlying 
event background. Figure [^a) shows an example of a 
seed cell and three of its first neighbors forming a four 
cell cluster. 

The second algorithm uses the track projection point 
onto the HBD to form a cluster around it. The pointing 
resolution of a track to HBD is ~3 mm at pr ~ 0.5 GeV/c 
which is much smaller than the size of a pad. The algo¬ 
rithm allows only up to three cells in a cluster, depend¬ 
ing on the track projection position within the cell. If 
the track projection points to the middle part of the cell, 
only that cell is used, but if it points to the edge of a 
cell one or two additional neighboring cells are summed 
up in the cluster [^. The same pattern of fired cells 
shown in Fig. [^a) would result in a three cell cluster in 
the projection-based algorithm as illustrated in Fig.j^b). 
The projection-based algorithm results in a more precise 
selection of the true hit, less fake hits and less pick up of 
charge from underlying event fluctuations. 

This is especially important in the most central colli¬ 
sions. On the other hand, the limited cluster size trun- 








cates the charge information, resulting in a somewhat re¬ 
duced efficiency and less power to discriminate between 
single and double hits. Therefore both algorithms are 
utilized in a complementary way, the stand alone pro¬ 
viding a higher efficiency and better single to double hit 
separation and the projection-based providing a better 
rejection of fake hits. 



FIG. 5. (Color online) (a) Stand-alone cluster formed by a 
seed cell (red) and three of its first neighbors resulting in 
a four cell cluster. Fired cells are colored, (b) The same 
pattern results in a three cell cluster with the projection-based 
algorithm that uses the projection point of an electron track 
onto the pad plane. 


E. Acceptance 


1. Acceptance during 2010 run 


As mentioned in Section IIB the PHENIX central arm 


magnets were operated in the H— configuration during 
the 2010 run. Compared to the standard -I—I- magnetic 
field configuration of PHENIX, the H— configuration has 
an increased acceptance for low pT tracks of about 20%. 

Charged particles are bent in the azimuthal direction, 
(p, by the magnetic field. Because the DC and RICH are 
needed to reconstruct the tracks and select the electron 
candidates, the azimuthal electron acceptance depends 
on their charge and px and on the radial location of each 
detector subsystem. We define the ideal track acceptance 
of the PHENIX detector in the H— field configuration by 
the following set of conditions: 


0min ^ 00 Q ^ 0max 

( 3 ) 

PT 


</>min < Po+q < (/>max 

( 4 ) 

PT 


^min ^ ^max 

( 5 ) 


and 0niax=l-92 rad are defined by the PHENIX central- 
arms pseudorapidity acceptance |? 7 | < 0.35. One of the 
arms covers the azimuthal range from (p^nin = ~ 

to (pmax = ^TT and the other from ^min = ig’’’ to 

'/'max = ffTr. The results shown in Section indicated 
as “in the PHENIX acceptance”, refer to the results fil¬ 
tered according to this parametrization of the ideal ac¬ 
ceptance. 


2. Fiducial cuts 

Several fiducial cuts are applied to remove inactive ar¬ 
eas of subsystems or areas with intermittent response, in 
order to homogenize the detector response over sizable 
fractions of the run time. Regarding the operation of the 
drift chamber, the entire 200 GeV Au-|-Au data set is di¬ 
vided into five groups, with fiducial cuts applied to each 
group separately such that inside each group the drift 
chamber has a stable active area. The nonactive DC ar¬ 
eas correspond to 19%-31% of the total DC acceptance, 
depending on the run group. 

Fiducial cuts are also applied to the HBD to exclude 
tracks pointing to one inactive module out of the 20 mod¬ 
ules of the HBD. Another fiducial cut removes conversion 
electrons originating from the HBD support structure, 
which are strongly localized in (p near the edges of the 
acceptance. Other fiducial cuts are applied to remove 
inactive or low efficiency areas in PCI and EMCal. 

In summary, the ideal PHENIX acceptance is reduced 
by the fiducial cuts by an amount that varies between 
32% and 42%, depending on the run group, with an av¬ 
erage of 36% for all selected runs. 


III. ANALYSIS 


This section describes the basic steps of the Au-|-Au 
data analysis. It is organized as follows. The data set 
and event selection cuts are presented in subsection |HI A| 
Subsection IHI Bl describes the track reconstruction. The 
methods applied to identify electrons are presented in 


detail in subsection |HIC| and the cuts applied to elec¬ 
tron pairs are explained in subsection mmi A detailed 
account of the various background sources and their sub¬ 
traction is provided in subsection [III E[ Next we present 
the raw spectra and corrections (subsection HI F) and 
discuss the systematic uncertainties (subsection [IIIG l. 
In the final subsection |HI H| we discuss a second indepen¬ 
dent analyses used as a cross check of the main analysis. 


A. Data set and event selection 


for tracks originating at z=0 with charge q, transverse 
momentum px and emission angles po and Oq. k^c = 
0.060 radxGeV/c and fcRicH = 0.118 radxGeV/c are the 
effective azimuthal bends to the DC and the RICH, re¬ 
spectively. The polar angle boundaries of 0niin=l-23 rad 


The Au-I-Au collision data at = 200 GeV were 

collected during 2010. Collisions were triggered using the 
beam-beam counters, with the MB trigger condition (see 
subsection HA). 
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The centrality is determined for each Au+Au collision 
from the sum of the measured charge in both BBCs com¬ 
bined with a Glauber model of the collision m as de¬ 
scribed in Ref. [H]. In this analysis, the data sample is 
divided into five centrality classes: 0%-10%, 10%-20%, 
20%-40%, 40%-60% and 60%-92%. The average num¬ 
ber of participants (A^part) and collisions {Ncoii} together 
with their systematic uncertainties associated with each 
centrality bin are summarized in Table [TTl 

TABLE II. Average values of the number of participants 
(A^part) and number of collisions (ALoii) for Au-|-Au collisions 
at = 200 GeV with the corresponding uncertainties. 

The values are derived from a Glauber calculation Eniiii]. 


Centrality 

(A^part )(sySt) 

(Mcoii) (syst) 

0 %-10% 

324.0 (5.7) 

951.1 (98.6) 

10 %-20% 

231.0 (7.3) 

590.1 (61.1) 

20%-40% 

135.6 (7.0) 

282.4 (28.4) 

40%-60% 

56.0 (5.3) 

82.6 (9.3) 

60%-92% 

12.5 (2.6) 

12.1 (3.1) 

0%-92% 

106.3 (5.0) 

251.1 (26.7) 


The data were recorded with an online vertex selec¬ 
tion of either ±20 cm (narrow vertex) or ±30 cm (wide 
vertex). The former selection was applied to the data 
recorded at the beginning of each store, when the lu¬ 
minosity was relatively high. For the latter selection, 
an additional-offline vertex cut of 30 < z < 25 cm was 
applied. This asymmetric cut is needed to avoid the in¬ 
creased yield of conversion electrons originating from the 
side panels of the HBD. These cuts resulted in 1.8 x 10® 
events with the narrow-vertex selection, 3.8 x 10® events 
with the wide-vertex selection, and a total of 5.6 x 10® 
MB events. 


B. Track reconstruction 

Charged particle tracks are reconstructed in the cen¬ 
tral arms using the DC and PCI [52]. The procedure 
assumes that all tracks originate from the collision ver¬ 
tex. Each reconstructed track is then projected onto the 
other detectors, RICH, EMCal, TOE and HBD, and the 
projection points are associated to reconstructed hits in 
these detectors. 

After a track is reconstructed, the initial momentum 
vector of the track at the z vertex is calculated. The 
transverse momentum pT is determined by measuring the 
angle a between the reconstructed particle trajectory and 
a line that connects the 0 -vertex point to the particle tra¬ 
jectory at a reference radius R = 220 cm. The angle a is 
approximately proportional to charge/pT- In the reverse 
field configuration used in the 2010 run, the momentum 
resolution is found to be 1.6% at pT = 0.5 GeV/c. 


C. Electron identification 

1. Detectors and variables used for electron identification 

For electron identification, the present analysis uses 
the HBD along with the central arm detectors RICH and 
EMCal and the time-of-flight information from the TOF- 
E detector and the EMCal. The relevant variables for 
electron identification from these detectors are: 

nO: number of hit PMTs in the RICH in the ex¬ 
pected range of a Cerenkov ring. 

disp: distance between a track projection and its 
associated ring center in the RICH. 

chi2/npe0: a x®-like shape variable of the RICH 
ring associated with the track per npeO, the number 
of photoelectrons measured in the ring. 

emcsdr: distance between the track projection 
point onto the EMCal and the associated EMCal 
cluster, measured in units of standard deviation of 
the momentum dependent matching distribution. 

prob: probability that the EMCal cluster is of elec¬ 
tromagnetic origin, based on the shower shape. 

dep: variable quantifying the energy-momentum 
matching for electrons. It is dehned as dep = 
^ , where E is the energy measured by the 
EMCal, p is the track momentum and aE/p is the 
momentum-dependent standard deviation of the 
Gaussian-like E/p distribution. 

stof(PbSc) and stof(TOF-E): time-of-flight de¬ 
viation from the one expected for electrons mea¬ 
sured by either the EMCal-PbSc or the TOF-E de¬ 
tector, converted in units of standard deviation of 
the Gaussian-like time-of-flight distribution. 

hbdcharge(P), hbdsize(P): cluster charge and 
size from the HBD projection-based algorithm. 

hbdid: reduced cluster charge threshold from the 
projection-based algorithm. This is the threshold 
of the hbdcharge(P) variable, that has been tuned 
to reduce the number of the nongenuine HBD hits 
by a hxed factor. E.g. by requiring hbdid>10, the 
number of the nongenuine HBD hits is reduced to 
l/IO of the initial number. These thresholds are 
tuned depending on event multiplicity and HBD 
cluster size. 

maxpadcharge(S): charge of the single pad with 
largest charge in the cluster of the stand-alone al¬ 
gorithm. 

hbdcharge(S), hbdsize(S): cluster charge and 
size from the stand-alone algorithm. 

First, electron candidates are selected from the total 
sample of tracks that contains mostly hadrons. This is 
accomplished by applying very loose cuts such as nO > 0, 
which requires at least one fired PMT around the track 
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projection in the RICH and E/p > 0.4 which rejects the 
tracks that strongly deviate from the expected E/p of 
~ 1. The sample of electron candidates selected in such a 
way comprises the signal electrons, background electrons 
(mostly conversions from the HBD back plane), and a 
relatively large number of misidentified hadrons. 


2. Exclusion of RICH photo-multipliers 

The RICH detector in PHENIX uses spherical mir¬ 
rors to project the Cerenkov light created by electrons 
in the radiator gas onto the PMT plane. As a conse¬ 
quence of this mirror geometry, parallel tracks after the 
field are projected to the same point in the PMT plane. 
In other words, if a hadron track is parallel to an electron 
track that produces a genuine response in the RICH, the 
hadron will appear to have the same response as the elec¬ 
tron and thus it will be misidentified as an electron. Fig¬ 
ure shows a typical example of this ring sharing effect. 
In this example, an electron-positron pair is generated 
by a photon conversion in the HBD backplane. After the 
magnetic field, a hadron track is parallel to the positron 
track. Consequently, the hadron and the positron share 
the same photomultipliers in the RICH detector and the 
hadron is misidentified as an electron. 



h- 


FIG. 6. (Color online) Illustration of a case leading to ring 
sharing in the RICH detector. The hadron track parallel to 
the positron track after the magnetic field will be misidentihed 
as an electron. 

This ring sharing effect occurs because the RICH re¬ 
construction algorithm allows multiple use of fired PMTs 
by different tracks. The ring sharing is a significant 
effect. In the 2010 run, the majority of electrons are 
generated by 7 conversion in the HBD backplane. Al¬ 


though these conversions can successfully be rejected by 
the HBD, their response in the RICH remains and there 
is some probability that the misidentified hadron will also 
remain in the pool of electron candidates. 

To reduce PMT sharing by different tracks in the 
RICH, the original RICH algorithm is modified. The 
PMTs fired by electrons that are clearly identified as 
background electrons, are removed, the ring reconstruc¬ 
tion algorithm is re-applied and new nO, npeO, disp, 
variables are derived. These background electrons are 
mainly conversion electrons from the HBD backplane, 
electron tracks pointing outside the HBD acceptance, 
electrons produced by conversion on the HBD support 
structure or low pp electrons with px < 200 MeV/c. 


3. The neural networks 

After the initial rejection of nonsignal electrons and 
the reduction of the ring sharing effect, the sample of 
electron candidates is still highly contaminated by back¬ 
ground electrons and misidentified hadrons. A standard 
procedure to increase the purity of the electron sample 
would be to apply a sequence of one-dimensional cuts 
on all or some of the fourteen variables listed above. 
However, such a procedure results in a large efficiency 
loss that becomes significant in the e“'"e“ pair analysis 
where the pair efficiency is approximately equal to the 
single track efficiency squared. In this analysis we im¬ 
plement instead a multivariate approach that is based on 
the neural network package TMultilayerPerceptron from 
ROOT [42] . 

The neural network comprises three layers: the input 
layer, the hidden layer and the output layer. The input 
layer is composed of all the input variables normalized 
to have their values between 0 and 1. The hidden layer 
comprises a selected number of neurons and the output 
layer comprises a single output variable. The number 
of neurons in the hidden layer determines the ability of 
the neural network to distinguish between the signal and 
the background, but this ability saturates with increasing 
number of neurons. For each neural network, we make 
sure that the number of neurons is sufficiently large to 
provide the best possible performance, typically 10-15 
neurons. In addition, we make sure that a sufficient num¬ 
ber of tracks is selected for the training sample, such that 
the performance of the neural network does not depend 
on the training statistics.The neural network output is a 
single probability-like variable, in which values closer to 
1 mostly correspond to signal, while values closer to 0 
mostly correspond to background (examples of the neu¬ 
ral network output distributions will be shown below). 
By selecting the tracks above a certain threshold, we can 
reject most of the background while keeping a large frac¬ 
tion of the signal. 

We use three different neural networks specially 
trained on subsets of the large list of elD variables to re¬ 
ject (i) hadrons misidentified as electrons in the central 
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arms (NNh), (ii) background electrons which are mostly 
HBD backplane conversions (NNe) and (ill) double hits 
in the HBD (NN^j). In this way we basically have three 
handles to separately treat each type of background. The 
neural networks learn to distinguish the signal and the 
background on well defined samples. The first two neural 
networks, NN^ and NNe, are trained on HUING events. 
The third neural network NNj is trained on a sample 
of single particle event simulations, (j) decays for 

single response and 7r° —>■ je^e~ Dalitz decays for dou¬ 
ble response. The training is done separately for each 
centrality bin in order to properly treat the multiplic¬ 
ity effects. For centralities > 40%, we use the neural 
network trained for the 20%-40% centrality bin, where 
the statistics of the training sample is higher. This is 
justified because already in the 20%-40% centrality bin, 
multiplicity effects are unimportant and the separation 
between signal and background is good. The training is 
also done separately for the three cases of time-of-flight 
information (TOF-E, PbSc-TOF, no time-of-flight infor¬ 
mation). 

The simulated events are passed through a GEANT sim¬ 
ulation of the PHENIX detector and through the same 
reconstruction code that is used for the data analysis. 
They are divided into two samples. One is used for train¬ 
ing purposes and the other one to monitor the neural net¬ 
work output. The simulated events are not used to de¬ 
termine absolute efficiencies (those are determined from 
simulation as discussed later in Section IIIF[ They are 
used only for training and monitoring purposes and the 
HUING events are particularly valuable in this respect. 
They allow us to assess the origin and relative magni¬ 
tude of the various background sources at each step of 
the electron identification chain, as well as the neural net¬ 
work performance in its ability to reject the background 
while preserving the signal. Details of the three neural 
networks are given below. 


4- Hadron rejection 

The first neural network, NNh, aims at reducing the 
hadron contamination. It exploits the information from 
all the relevant detectors, HBD, RICH, EMCal and TOF- 
E. The signal (S) for the training of NNh comprises elec¬ 
tron tracks originating at the collision vertex, whereas 
the background (B) comprises all the remaining misiden- 
tified hadron tracks in the sample. 

Eigurej^shows the output values of NNh for the hijing 
monitoring sample (red line) and also shows the output 
of NNh applied on real data (black line). The truth in¬ 
formation from the hijing events in terms of signal and 
background is shown separately. It should be noted that 
in the hijing monitoring sample, all electron tracks are 
considered. The signal comprises the genuine electrons 
excluding the HBD backplane conversions and the back¬ 
ground is all remaining tracks. 



NN^ output 


FIG. 7. (Color online) Comparison of the output values of the 
neural network NNh for the 0%-10% centrality bin applied 
to the HIJING monitoring sample (red line) and to real data 
(black line). The figure also shows the signal (green) and 
the background (blue) components of the hijing simulation. 
The arrow represents the average final cut selected by the cut 
optimization procedure. See text in Section [ill C 7| 


5. Background electron rejection 

After rejecting hadrons in the previous step, the dom¬ 
inant background in the electron sample comes from the 
conversions in the HBD backplane that were not rejected 
by the conservative process described in |HIC2| Because 
these conversions do not leave a signal in the HBD they 
can be recognized and rejected if the tracks do not have 
a matching HBD response. The rejection capability is 
however limited by fluctuations remaining after the un¬ 
derlying event subtraction in the HBD. To provide the 
optimal rejection of the remaining backplane conversions 
we use a neural network, NNe, which is based on the 
HBD information reconstructed by both the stand-alone 
and the projection-based algorithms. The signal tracks 
for the training of NNe comprise all signal electrons re¬ 
maining after the previous step, while the background 
sample includes only the electrons originating from the 
HBD backplane. 

Figurej^shows the distribution of output values of NNe 
applied to the hijing monitoring sample (red line) and 
to data (black line). The signal and background compo¬ 
nents of the hijing simulation are shown separately. 


6. Double-hit rejection in the HBD 

After removing hadrons and backplane conversions as 
much as possible, the major sources of background are 
the beam-pipe and radiator conversions and electrons 
from Dalitz decays where only one track is recon¬ 
structed in the central arms. These electrons have a 
zero or very small opening angle and most of them lead 
to a double hit in the HBD. Double hits can be recog¬ 
nized using the HBD response reconstructed in parallel 
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FIG. 8. (Color online) Comparison of the output values of the 
neural network NNe for the 0%-10% centrality bin applied 
to the HIJING monitoring sample (red line) and to real data 
(black line). The figure also shows the signal (green) and 
the background (blue) components of the hijing simulation. 
The arrow represents the average final cut selected by the cut 
optimization procedure. See text in Section [ill C 7| 


double response by electrons from tt^ —> ye+e” Dalitz 
decays. The simulations are embedded into real HBD 
background events in order to take into account centrality 
dependent occupancy effects. 


7. Cut optimization 

The final selection of cuts on each neural network 
output variable is optimized using hijing events. The 
thresholds are varied separately to maximize the effec¬ 
tive signal, Because the statistics of the hijing 

samples are by far insufficient for a pair analysis, for the 
signal S we use the number of single electrons from charm 
decay per event, which is an easily identified signal in hi¬ 
jing, and for the background B we use the total number 
of electrons per event. The cut optimization is done sep¬ 
arately for each centrality class, for two ranges {pt 
< 300 MeV/c and pr > 300 MeV/c), for each cluster 
size, and for each TOF conhguration. The effective sig¬ 
nal for each setup is maximized subject to the following 
conditions: 


by both the stand-alone and the projection-based algo¬ 
rithms. The response is coupled in a neural network, 
NNd separately optimized for different BBD cluster sizes 
as well as centrality classes. The NNd cut is an implicit 
small opening angle cut given by the maximum cluster 
size which is of the order of 75 mrad. 


• The three types of TOF configuration (with PbSc 
timing information, with TOF-east timing infor¬ 
mation and without any timing information), have 
similar efficiencies with differences of less than 15%. 

• Badron contamination less than 5% for TOF-E and 
PbSc-TOF and less than 10% for the no-TOF case. 



FIG. 9. (Color online) The output of the neural network 
NNd for the recognition of single and double hits in the HBD. 
Single response (solid line) is provided by electrons from sim¬ 
ulated <j} —decays and double response (dashed line) 
by electrons from tt® —i 7 e''"e“ Dalitz decays. This example 
is for 30%-40% centrality and for a three cell cluster size. 
The arrow represents the average final cut selected by the cut 
optimization procedure. See text in Section |ra^ 


Figure shows the distribution of the output variable 
of the neural network NNd for the separation of single and 
double hits in the BBD. The single response is provided 
by electrons from simulated (j) —decays and the 


The arrows in Figs. [7][^ represent the average final cuts 
selected by the cut optimization procedure for these par¬ 
ticular cases. The final cuts produce an electron sample 
with small hadron contamination, of less than 5%, for 
all centralities. Strong cuts on the BBD are needed to 
achieve this small hadron contamination, resulting in a 
single electron efficiency of 25%-40% d epend ing on cen¬ 
trality, at Pt > 0.5 GeV/c (See Section BIF). 


D. Pair cuts 

The track selection criteria described above provide an 
electron sample with high purity. Bowever, besides these 
criteria which are applied on a track-by-track basis, this 
analysis implements a series of dielectron cuts, based on 
the pair properties. These cuts are needed in order to 
remove ghost pairs i.e. pairs correlated by the close prox¬ 
imity of tracks in one of the detectors. Such correlations 
cannot be described by the mixed background, by def¬ 
inition, therefore this part of the phase-space must be 
removed from both the foreground and the mixed back¬ 
ground. In the present analysis we remove the whole 
event, if such a pair is found, as was done in Ref. [23]. 
This procedure removes only ^2% more of the total pair 
yield than discarding the pairs, because the average pair 
multiplicity is relatively low. 
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The most prominent detector correlation comes from 
the ring sharing effect in the RICH detector, discussed in 
Section pll C 2[ which arises when two tracks are parallel 
after the magnetic field, with at least one of them being 
an electron. 

As mentioned above, the detector-correlated pairs are 
identihed by applying a cut on the physical proximity 
of the tracks forming a pair in every detector and the 
cut value is determined by the corresponding double hit 
resolution. In the RICH detector, the cut selects pairs 
whose rings are closer than 36 cm, which is twice the 
diameter of the RICH ring (~16.8 cm). In the EMCal, 
the cut removes a region of 2.5 x 2.5 towers around the 
hit. In PCI the pairs are selected for removal if their 
tracks are within 5 cm in z or 0.02 rad in (f). 

The effect of these three pair cuts on the like-sign and 
unlike-sign mass spectra is shown in Fig. 10 The like- 
sign yield close to ruee ^ 0 GeV/c^ is affected by all 
cuts. On the other hand, in the unlike-sign foreground 
spectrum, the cuts affect well localized regions producing 
two clearly visible dips. The dip at rUee ^ 0.25 GeV/c^ 
is created by the RICH pair cut and the dip at niee ~ 
0.15 GeV/c^ is created by the PGl pair cut. The EMGal 
pair cut removes yield around 0.20 GeV/c^, but the effect 
is small compared to the other two cuts. 

In addition to the RICH, EMCal and DC/PCl ghost 
cuts, a 100 mrad opening angle cut is applied to remove 
ghost pairs in the HBD. This is a proximity cut that 
translates to a distance of two cells in the pad readout 
and roughly corresponds to the double hit separation of 
the HBD. This cut affects the yield at niee 0 GeV/c^ 
in both the like-sign and unlike-sign mass spectra. 



0 0.1 0.2 0.3 0.4 0.5 0.6 


(GeV/c^) 

FIG. 10. (Color online) (a) Like-sign and (b) unlike-sign 
foreground spectra without any pair cuts (Black) and with 
RICH, EMCal and PCI pair proximity cuts (Blue) for MB 
events. 


E. Background Pair Subtraction 

Because the origin of the electron track candidates is 
not known, all electrons and positrons in the same event 

are paired to form the unlike-sign (EG-|_) and like-sign 

(EG_|_+ and FG _) foreground mass spectra. This gives 

rise to a large combinatorial background that increases 
quadratically with the event multiplicity. In addition to 
that, there are several background sources of correlated 
pairs. The evaluation and subtraction of the background 
is the crucial step in the analysis of dileptons in particular 
in situations, like the present one, where the S'/H is at the 
sub-percent level. In this section, we describe in detail 
the various sources contributing to the background and 
the methodology used to evaluate each of them. 


1. Background sources 

The unlike-sign foreground spectrum FG^ _contains, 

in addition to the physical signal (S'), a large background 
comprising the following sources: 

• Uncorrelated combinatorial background (GB): It 
arises from the random combinations of electrons 


and positrons originating from different parent par¬ 
ticles and is an inherent consequence of pairing all 
electrons with all positrons in the same event. The 
combinatorial background accounts for most of the 
total background, more than 99% in the most cen¬ 
tral collisions and more than 90% in peripheral 
collisions. The two electron tracks of combinato¬ 
rial pairs are uncorrelated. However, they carry 
a global modulation induced by the collective flow 
of each individual collision. The evaluation of the 
combinatorial background together with the flow 
modulation is described in detail in the following 
subsection. (See Section HIE2 ) 


• Correlated background pairs. There are three dif¬ 
ferent sources of correlated background pairs: 

— Gross pairs {CP): A cross pair can be pro¬ 
duced when there are two e+e“ pairs in the 
final state of a single meson decay. One such 
case is 7 r° —>■ e+e “7 —>■ e+e^e+e”. The pair 
formed by an electron directly from 7 r° and 
a positron from 7 conversion does not come 
from the same parent particle but it is a corre¬ 
lated pair through the same primary particle. 
(See Section HIE3 ) 
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— Jet pairs {JP)'- The jet pairs are produced by 
two electrons generated in the same jet or in 
back-to-back jets. (See Section HIE4 ) 

— Electron-hadron pairs {EH)-. Whereas the 
previous two sources of correlated pairs are of 
physics origin, the electron-hadron pairs are 
an artifact that results from residual detec¬ 
tor correlations that cannot be handled by the 
pair cuts. (See Section HIE5 ) 


One can then write: 

EG+_ = S + CB+_ + CP+_ + JP+_ + EH+_ (6) 

All the background sources listed above form the yield 
of the like-sign foreground mass spectra FG+.\. and 
EG _There is no signal in these spectra with the ex¬ 
ception of a very small contribution of and e“e“ 

pairs from bb decays {BB). So one can write: 


FGj^j^ — GBj^j^ CPj^j^ -|- J T++ -f -|- (7) 


EG— = GB__ + GP— -k JP— + EH__ + (8) 

Usually the like-sign pairs are subtracted from the 
unlike-sign pairs to obtain the signal. This is a convenient 
approach in a detector with 2tt azimuthal coverage, which 
ensures that the uncorrelated background is charge sym¬ 
metric, under the assumption that the correlated back¬ 
ground is also charge symmetric, i.e. it produces the 
same yield and mass distribution of like and unlike pairs. 
These conditions are not met in the present situation. 
The two central arm configuration of the PHENIX de¬ 
tector results in a substantial acceptance difference be¬ 
tween like and unlike-sign pairs. Furthermore, the like- 
sign pairs contain a small signal component from bb de¬ 
cays that needs to be calculated separately. Finally, as 
shown below, the electron-hadron pairs are not charge 
symmetric. For these reasons, in this analysis we adopt 
a different approach in which each source is evaluated 
separately for a quantitative understanding of the like- 
sign yield. Once this is demonstrated, the background 
sources, CB,GP^JP and EH are subtracted from the 
inclusive foreground unlike-sign spectrum in order to ob¬ 
tain the mass spectrum of the signal pairs. The following 
subsections outline the evaluation of the various back¬ 
ground sources. 

The BB contribution which is part of the signal 
is needed only for the quantitative evaluation of the 
like-sign spectra. The contribution is calculated using 
MC@NLO (See Section IV for details), which generates 


both like-sign and unlike-sign contributions from BB. 
The small like-sign contribution from DD is neglected. 


events but with similar characteristics are combined into 
pairs. In this analysis, all events are classified into 11 
bins in 2 vertex between —30 cm and -1-25 cm, and 10 
bins in centrality between 0% and 92%. 

In principle, the event mixing technique is expected 
to reproduce the shape of the combinatorial background 
with great statistical accuracy, because one can mix as 
many events as needed to reduce the statistical uncer¬ 
tainty to a negligible level. In fact it does not reproduce 
the shape. There is a small difference between the fore¬ 
ground combinatorial background and the mixed event 
background. The former is affected by the elliptic flow 
which is intrinsic to heavy ion collisions, whereas the lat¬ 
ter is obtained by randomly picking up two tracks from 
different events and thus on the average does not have 
any flow effect. 

To take into account the effect of flow in the mixed- 
events, one could make reaction plane bins, in addition 
to the vertex and centrality bins, so that only events with 
similar reaction plane are mixed. However, the method is 
limited by the reaction plane resolution and in PHENIX, 
the latter is not sufficient to reproduce the shape of the 
foreground combinatorial background. Instead, in the 
present analysis, a weighting method, based on an an¬ 
alytical calculation of the flow modulation, is used to 
account for the flow effects in the mixed events. 

If particles are generated according to the following 
distribution function: 


1 J-2z; 2 cos 2((/) — ■!/;), (9) 


where (j) is the particle emission angle in azimuth, ip is 
the reaction plane angle and V 2 is the elliptic flow coeffi¬ 
cient, then random pairs formed from these particles are 
distributed as (See Appendixfor the derivation): 

P{(pa -(t)b) = l + 2v2,aV2,bCOS2{cl)a “ (pb), (10) 


where 4>a{b) is the azimuthal emission angle and V2,a{b) 
the elliptic flow of the two particles forming the pair. 

In the weighting method, each mixed background pair 
is weighted by Eq. (10). The V 2 values of inclusive elec¬ 
trons are determined from the present data prior to the 
pair analysis as a function of centrality and electron pT 
using the reaction plane method [43]. Exactly the same 
cuts as in the data analysis are used in the V 2 calculation. 
The obtained V 2 values are in very good agreement with 
the inclusive electron V 2 values reported in Ref. [44j . 

We use a Monte-Carlo (MC) simulation to evaluate 
the method. The simulation generates electrons and 
positrons following a Poisson distribution with a mean 
value of three The particles are uniformly distributed 


2. Combinatorial background (CB) 

The combinatorial background is determined using the 
event mixing technique, in which tracks from different 


^ There is not much meaning to the mean value of 3 of the Poisson 
distribution. It is a convenient choice to have one pair per event 
with a high probability. 
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in pseudorapidity between ±0.35 and their momentum 
distribution is taken from data. The azimuthal emis¬ 
sion angle (j) is determined according to the distribution 
1 ± 2 u 2 COs 2 (^ — '0), where 0 is the reaction plane an¬ 
gle, which is uniformly distributed between ±^. The V 2 
values are taken from the 20%-40% centrality bin. The 
tracks that pass the PHENIX acceptance filter are used 
in the pair analysis. 
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FIG. 11. (Color online) Foregronnd to mixed background ra¬ 
tio of (a) like-sign and (b) unlike-sign mass spectra ratio in 
a MC simulation. The foregronnd is generated with flow, 
whereas the mixed events are produced without flow i.e. us¬ 
ing a simple mixed-event technique (squares) and with flow 
modulation using the weighting method (circles). 


Figure [m shows the ratio of the foreground to mixed 
background mass spectra. The squares correspond to the 
simple mixed-event technique without correcting for flow. 
We can see that in this approach the ratio is not flat, i.e. 
the foreground shape is not reproduced by the mixed 
background shape. The circles correspond to the weight¬ 
ing method. The ratio is completely flat over the en¬ 
tire mass range demonstrating that the weighting method 
properly accounts for the flow modulation. 


3. Cross pairs (CP) 

Cross pairs can be produced when a hadron decay pro¬ 
duces two e“''e“ pairs in the final state. The following 
hadron decays and subsequent photon conversions lead 
to cross pairs: 


TT^ ^ ef Cl ^ ^ ef 63 ( 11 ) 

7r° ^ 7172 ^ eC eje0 (12) 

7 )—i eO j—i 62 (13) 

77 ^ 7172-)> 6)^6]" eje0 (14) 


The cross combinations give rise to two unlike-sign 
pairs (e ^62 and e^e)”) as well as two like-sign pairs 
(e^e^ and e)”e 0 ) that are not purely combinatorial, but 
correlated via the or 77 mass and momentum. There¬ 
fore, this contribution is not reproduced by the event¬ 
mixing technique. 

To calculate the cross pairs, we use EXODUS (see Sec¬ 
tion IV I to generate 7r° and 77 with the following input 
parameters: 


• Flat-vertex distribution within |z| < 30 cm. The 
hnal results are weighted to restore the measured 
vertex distribution. 


• Flat pseudorapidity distribution within I 77 I < 0.6 
and uniform in 0 within 0 < 0 < 27r. 


• Momentum distributions based on PHENIX mea¬ 
surements (see Section IV). 


The generated 7r° and 77 are passed through a GEANT 
simulation of the PHENIX detector. By selecting recon¬ 
structed cross pairs, one can determine the shape of the 
cross-pair invariant mass spectrum. The spectra are then 
absolutely normalized using the rapidity density values 
dN^^o/dy and dN^^/dy as a function of centrality, sum¬ 
marized in Section |IV| The absolutely normalized mass 
spectra of cross pairs for the 0 %- 10 % centrality bin are 
shown in Fig. [T^ 


4- Jet pairs (JP) 

The jet pairs are produced using the pythia 6.319 code 
with CTEQ 5 L parton distribution functions [20]. The fol¬ 
lowing hard quantum-chromodynamics (QCD) processes 
are activated [50] : 


A similar MC study was performed to evaluate whether 
triangular flow 773 also induces shape distortion of the 
mass spectrum. For the most central collisions, where 
773 is comparable to V2 at high pp [IS], the simulations 
show that the 773 effect is at least one order of magnitude 
smaller than for V 2 and we thus ignore triangular flow 
in the determination of the combinatorial background 
shape. 


. MSUB II: /,/, ^ hf, 

. MSUB 12: /J, ^ fj, 
. MSUB 13: fj, ^ gg 
. MSUB 28: Ug ^ Ug 
. MSUB 53: gg ^ h~fk 
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FIG. 12. Absolutely normalized (a) like-sign and (b) unlike-sign spectra of cross pairs {CP) from EXODUS and GEANT simulations 
for the 0%-10% centrality bin. The and r/ contribntions are shown separately. 


MSUB 68 : gg gg 


where g denotes a gluon, fij^k are fermions with flavor 
i, j, k and fi^j^k are the corresponding antiparticles. A 
Gaussian width of 1.5 GeV/c for the primordial kr dis¬ 
tribution (MSTP(91)=1, PARP(91)=1.5) and 1.0 for the 
K-factor (MSTP(33)=1, PARP(31)=1.0) are used. The 
minimum parton pr is set to 2 GeV/c (CKIN(3)=2.0). 
The z coordinate of the vertex position is produced uni¬ 
formly between ±30 cm and then weighted to reproduce 
the measured distribution. From the pythia output, 7r° 
and r] are extracted and passed through the GEANT simu¬ 
lator of PHENIX in order to generate the inclusive e“'"e“ 
pairs. 

In addition to the jet pairs we are interested in, the 
foreground pairs from pythia events contain also “phys¬ 
ical” pairs, cross pairs and combinatorial pairs. The 
“physical” pairs and cross pairs are excluded from the 
foreground pairs by requiring that the two electrons or 
positrons of the pair do not share the same particle in 
their history. The combinatorial background is statisti¬ 
cally subtracted using the event-mixing technique. The 
mixed event like-sign pairs are normalized to the fore¬ 
ground like-sign pairs in the range ~ t/2 , where 

is the difference in the azimuthal angle of the 


primary particles, or rj. Figure 
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shows the A(/)g”*"* 


distributions of the foreground pairs and the normalized 
mixed-event pairs. The excess yield around ~ 0 

represents the dileptons from the same jet whereas the 
excess yield at Ai^g”™ ~ tt corresponds to the dileptons 
from opposite or back-to-back jets. 

After subtracting the combinatorial background, the 
PYTHIA spectra are scaled to give the pion yield per p+p 
MB event . The scaling factor is determined such that the 
7r° yield in the pythia simulation matches the measured 
7 r° yield in p+p collisions [47j and found to be 1/3.9. 

The spectra need to be further scaled to obtain the 
jet contribution in Au±Au collisions for each centrality 
bin. This scaling is done following Ref. [IS]: an ee jet 
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FIG. 13. (Golor online) A())g’’™ (difference in the az¬ 
imuthal angle of the primary particles, 7r° or g) distributions 
of foreground and normalized mixed-event background like- 
sign pairs as obtained from the pythia simulations. 


pair originating from primary particles with momenta 
Pt,i and pt ^2 is scaled by the average number of binary 
collisions (Ncow) for each centrality bin, times RAAiPTp), 
times Iaa(j>t, 2 )- The same jet or opposite jet Iaa{pt, 2 ) 
values are applied depending on the pair opening angle. 
The absolutely normalized jet pair spectra for the 0%- 


10% centrality bin are shown in Fig. 14 


5. Electron-hadron pairs (EH) 

Even after applying the pair cuts described in Sec¬ 
tion |III D| electron-hadron pairs correlated through de¬ 
tector effects remain in the foreground pairs. An example 
of such an electron-hadron pair can be illustrated with 
the sketch of Figure discussed in Section El In 
this example, if both the positron and the mis-identified 
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FIG. 14. Absolutely normalized (a) like-sign and (b) unlike-sign spectra of jet pairs {JP) simulated by pythia and GEANT for 
the 0%-10% centrality bin. The near-side and away-side contributions are shown separately. 


hadron are detected, the pair is identified as a RICH 
ghost pair and the entire event is rejected by the RICH 
ghost pair cut as described in Section [Til D[ However, if 
the positron is not detected due to detector dead areas or 
reconstruction inefficiency, the pair formed by the elec¬ 
tron and the mis-identified hadron is not rejected and 
remains in the sample. This pair is not a combinatorial 
pair but correlated through the positron. Although the 
mis-identification of hadrons via hit sharing occurs in all 
detectors, the RICH detector is the dominant contribu¬ 
tor to these electron-hadron pairs. Therefore, only the 
RICH detector is considered as the source of such corre¬ 
lated pairs. 

We simulate electron-hadron pairs using electrons from 
TT^ and 7 ] simulations and hadrons from real events. The 
TT^ and r] simulations are the same ones that are used for 
the cross pair simulation. The hadrons from real events 
are all the reconstructed tracks that fail the elD cuts. 

The simulation is performed in the following way: 
First, a combined event is formed using electrons from 
one Dalitz decay of 7 r° or 77 generated with EXODUS and 
hadrons from a real event. Second, the information from 
their associated fired PMTs is merged and new rings are 
reconstructed. Using the new RICH ring variables, the 
regular analysis procedure, including elD cuts and pair 
cuts, is performed on the combined event. Finally, the 
pairs formed by the combination of an electron track 
from simulation and a hadron track from data are ex¬ 
tracted. The spectra are absolutely norm alized using the 
7 r° dN/dy values shown in Section IV The absolutely 
normalized electron-hadron pair spectra for the 0 %- 10 % 
centrality bin are shown in Fig. |15| Contrary to the cross 
pairs and the jet pairs where the like- and unlike-sign 
spectra have a very similar shape, the electron-hadron 
pairs exhibit a sizable difference between the like- and 
unlike-sign spectra. The yield of electron-hadron pairs 
has a strong centrality dependence. It increases by a 
factor of ^50 from peripheral to central collisions with 
respect to the tt^ rapidity density. This increase is mainly 


due to the expected scaling of the electron-hadron pairs 
with the square of the event multiplicity. 


6. Background normalization 


The cross pairs, jet pairs, electron-hadron pairs and 
bb decay pairs are absolutely normalized. The mixed 
event technique provides only the shape of the combi¬ 
natorial background. It needs to be normalized in order 
to be able to subtract the background and extract the 
signal. The only free parameters of the entire procedure 
are thus the normalization factors of the mixed event 

background like-sign spectra n/++ and nf _They are 

determined by normalizing the mixed event background 
yield ) to the foreground yield (iVFG++(__)), 

integrated over a selected region of phase space, after 
subtracting the correlated pairs integrated over the same 
region: 


nf++ 
nf — 


Nfg++ — Ncp++ — Njp++ — Neh++ — Nbb++ 
Nmix+^ 

NfG-- — Ncp__ — Njp__ — Neh-- — Nbb-- 
Nmix-- 


where Ncp^^^__^, Neh++^__) and 

_) are the integral yields of each source in 

the normalization region. The normalization region is 
a window in the azimuthal angular distance of the two 
tracks A^q- It needs to satisfy two competing condi¬ 
tions. On the one hand, a small normalization window 
containing only combinatorial pairs is preferred to avoid 
being affected by any residual yield (and systematic un¬ 
certainties) from the correlated background sources. On 
the other hand, a wide normalization window is required 
to reduce statistical uncertainty. The normalization 
windows used in this analysis for each centrality bin 
are shown in Table in together with the corresponding 
number of like-sign pairs {N^s = Nfg++ + Nfg _)■ 
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FIG. 15. Absolutely normalized (a) like-sign and (b) unlike-sign spectra of simulated electron-hadron pairs (EH) for the 
0%-10% centrality bin. See text for details. 


The region of small opening angles that correspond to 
small masses where the correlated pairs CP, JP and 
EH mostly contribute, is excluded in all centrality bins. 


nf _as: 

nf+- = a/ nf++ ■ nf— (19) 


TABLE III. Normalization window for each centrality bin. 
The number of like-sign pairs Nls in the window is also 
shown. 

Centrality 

Normalization window 

A(j>o 

Nls 

0 %-10% 

0.7 - 3.14 

5.1M 

10 %-20% 

0.7 - 2.1 

I.IM 

20%-40% 

0.7 - 2.1 

660K 

40%-60% 

0.9 - 2.1 

48K 

60%-92% 

0.9 - 2.1 

3K 


The combinatorial background in Eqs. 0 and Q is 
thus given by the normalized mixed-event background: 


CP++(mee) = nf++ ■ MIX++{mee) (15) 

CB {niee) = nf • MIX (Wee) (16) 

As long as electrons and positrons are produced in 
pairs and these pairs are uncorrelated, the total unlike- 
sign combinatorial background yield is the geometric 
mean of the total like-sign combinatorial yield, indepen¬ 
dent of single electron efficiency and acceptance [23]: 

CP+_ = 2^CB++ ■ CB__ (17) 

A similar relation holds true for the integral yields of the 
mixed-event background: 

MIX+_ = 2^MIX++ ■ MIX__ (18) 

The normalization factor n/-|_of the unlike-sign mixed 

event background is thus deduced from the normaliza¬ 
tion factors of the like-sign mixed background, n/_|__|_ and 


In the present analysis, the square root relation, Eq. 
0 . is violated by two independent factors. First, the 
relation does not hold true when pair cuts are applied to 
the spectra because pair cuts affect differently the unlike- 
sign and like-sign spectra. Second, elliptic flow induces 
an inherent distortion of the square root relation. Flow 
does not create or destroy particles. It only affects their 
azimuthal distribution and therefore in a perfect 27r de¬ 
tector there is no effect and Eq. 0 is obeyed. How¬ 
ever, in the case of the PHENIX detector, which is not 
a 27r detector, the relation is violated as demonstrated 
in Appendix Relation (191 can still be used provided 
that the violation is the same in the data and the mixed 
events. In the present analysis, we make sure that this is 
the case. We start from a situation in which the mixed 
events satisfy Eq. (18). We then apply to the mixed 
events the pair cuts, exactly as to the foreground events, 
and the flow modulation using a weighting factor pro¬ 
cedure that is based on an exact analytical calculation. 
Thus we make sure that Eq. (191 is still valid. 


7. Quantitative understanding of the background 

To illustrate our understanding of the background in 
quantitative terms. Fig. shows a comparison of the 
MB mass spectra for the foreground and the calculated 
background like-sign pairs. 

The top panel shows the foreground like-sign mass 
spectrum (open circles) together with the various back¬ 
ground components discussed above (the normalized 
combinatorial background, and the absolutely calculated 
cross pairs, jet pairs and e-h pairs) and the bb pairs calcu¬ 
lated as described in SectionjlVj The bottom panel shows 
the ratio of the foreground like-sign spectrum to the sum 
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FIG. 16. (Color online) (a) Measured like-sign spectrum 
(open circles) together with the calculated background com¬ 
ponents (histograms) for MB events, (b) Ratio of the like-sign 
spectrum to the sum of all the background components. 


of all the background components. Similar comparisons 
for the five centrality bins used in this analysis are shown 
in Fig. 

In general the background is well reproduced both in 
shape and magnitude. In particular, for the most cen¬ 
tral bins, the background is reproduced with sub-percent 
accuracy. There are, however, a couple of regions where 
the ratio foreground/background is different from one. 
There is a deviation of the order of a few percent at 
masses ruee < 100 MeV/c^. This is clearly visible in the 
three most central bins. A number of factors could be 
responsible for this deviation, such as scale errors in the 
cross pairs or the jet pairs. However, in this mass region 
the signal to background ratio is relatively good as shown 
in Fig. and a deviation of the order of a few percent 
in the background is negligible. There also seems to be 
a deviation at mee > 1 GeV/c^ for the 10%-20% and 
20%-40% centrality bins. This deviation could indicate 
underestimations of the flow or the back-to-back jet con¬ 
tributions, due to the precision in these measurements, 
or the existence of an additional correlation that is not 
taken into account in any of the calculated background 
components. To be conservative, this deviation is con¬ 
sidered as evidence of unsubtracted background and its 
magnitude is assigned as a mass dependent systematic 
uncertainty of the signal. 

Figure[T^shows the MB mass spectra of the foreground 
unlike sign events (FG^_), the calculated total back¬ 
ground (BG^_) and the raw signal obtained by their 



FIG. 17. Ratios of the like-sign foreground spectrum to the 
sum of all the background components for the five centrality 
bins used in this analysis. 


subtraction. The signal to background ratio is shown in 
the bottom panel. This result will be discussed in refer¬ 
ence to previously published PHENIX results in Section 
IVCT] 



FIG. 18. (a) MB mass spectra of the unlike sign foreground 

events (FG^_), the calculated total background (BG^_) and 

the raw signal S. (b) The signal to background ratio. 
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F. Raw Spectra and Efficiency Corrections 


Figure [I^ shows the raw mass spectra, obtained after 
subtracting the pair background, for the five centrality 
bins of this analysis. 

> 

0 

CD 

> 

CD 

2 
0 


10 ‘® 


10 -^^ 
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FIG. 19. (Color online) Raw mass spectra for the five cen¬ 
trality bins. 


To obtain the invariant mass spectrum inside the ideal 
PHENIX acceptance, the e“''e“ raw mass yield is cor¬ 
rected for reconstruction efficiency effects according to: 

dN _ 1 N{mee) 1 ^ 20 ) 

dmee ^evt ArUee 

where Ng^t is the number of events, N{mee) is the num¬ 
ber of e“'"e“ pairs with invariant mass Wee and Arn-ee is 
the mass bin width. is the total pair reconstruction 

efficiency that includes the elD efficiency of the neural 
networks, losses incurred by dead or inactive areas in the 
detector, pair cut losses and detector occupancy effects. 
The total pair reconstruction efficiency can thus be 
written as: 


total elD ^ live ^ ghost ^ mult 

^pair ^pair ’ ^pair ’ ^pair ’ ^pair 


( 21 ) 


where ep{f^ is the e+e” pair reconstruction efficiency in¬ 
cluding the efficiency of ail the electron identification cuts 
and the HBD doubfe-hit rejection cut, Cph® is the pair ef¬ 
ficiency from the detector active area with respect to the 
ideal PHENIX detector acceptance, reflects the ef- 

hciency loss due to the pair cuts that remove ghost pairs 


in the various detectors (see Section IIID|) and is 


the multiplicity dependent efficiency loss discussed below 
in this subsection. 

The single electron reconstructio n effi ciency, defined as 

vs pt for the five 


-el D ^mult 


pair 


shown in Fig. 
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centrality bins. This efficiency is not actually used in the 
analysis. It is shown here for illustration purposes. The 
change of efficiency below 0.3 GeV/c arises from the cut 
optimization in two px ranges (see Section IIIC7). 



FIG. 20. (Color online) Single electron reconstruction effi¬ 
ciency vs. Pt for the five centrality bins. 

The product ep^f^ • ■ ep!j°/* is determined as follows. 

A cocktail of all the known hadronic sources contribut¬ 
ing to the e+e” pair spectrum is generated within I 77 I < 
0.6 and 27r in azimuthal angle. Details about the vari¬ 
ous sources of the cocktail are given in Section |IV] The 
cocktail is passed through a full GEANT simulation of the 
PHENIX detector and analyzed in the same way as 
the data, including elD cuts, fiducial cuts and pair cuts. 
The resulting output is referred to as the reconstructed 
cocktail. The ratio of this reconstructed cocktail to the 
generated cocktail filtered through the ideal PHENIX ac¬ 
ceptance (but without momentum smearing), gives the 
product This correction is derived in 

the two dimensional space of mass-pair px- 

Special care is taken to tune the simulations to the data 
to ensure that the detector response in the simulations is 
the same as in real data for all the subsystems involved 
in the analysis. As an example, Eig. shows a compar¬ 
ison of a few electron identihcation variables in data and 
simulations. For this comparison we use a clean sample 
of electrons provided by fully reconstructed tt® Dalitz de¬ 
cays with an opening angle larger than 100 mrad from 
the 60%-92% centrality bin where the occupancy effects 
are very small and can be ignored. The elD variables of 
the two tracks from these pairs are compared to those of 
TT^ —>■ e+e“ 7 simulations. 

The HBD occupancy effects are taken into account by 
embedding the HBD hits from the cocktail simulation 
into real HBD events, and thus are included in the prod- 
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TABLE IV. Efficiency loss due to detector occupancy in the 
central arms and to the tagging of RICH PMTs dis¬ 

cussed in Section C 2| for the five centrality bins used in 
this analysis. 



O 

1 — 1 

1 

o 

10%-20% 

Centrality 

20%-40% 

o 

1 

o 
^ ' 

60%-92% 

embed 

^pair 

TPMT 

^pair 

0.53 

0.88 

0.65 

0.92 

0.76 

0.94 

0.86 

0.98 

0.95 

1.00 


Figure shows the total pair reconstruction efficiency 
^pair* for pair px within 0.8-1.0 GeV/c for each centrality 
bin. 


FIG. 21. (Color online) Comparison of electron identifica¬ 
tion variables in data (black) and in simulations (red). The 
variables are described in Section Fill Cl electrons in data and 
simulations are from fully reconstructed tt® Dalitz decays with 
opening angle larger than 100 mrad. 


net • ep™if. • ep!j“*. There are two other occupancy 
effects in the central arms that need to be taken into 


account and are included in Eq. (21) by the additional 


multiplicative factor The first one is the decrease 

of track reconstruction efficiency as the detector occu¬ 
pancy increases with centrality. This loss is referred to 


as and is determined by an embedding procedure. 

Electrons from </> decays that are reconstructed in sin¬ 
gle particle simulations, are embedded into real Au-|-Au 
events. Then the embedded events are run through the 
full reconstruction software chain and analyzed in exactly 
the same way as the data. The embedding efficiency for 
single tracks is determined as the ratio of the num¬ 

ber of reconstructed electron tracks from embedded data 
to the number of embedded tracks. The pair embedding 
efficiency is calculated as the square of the single track 


embedding efficiency, e® 


= (< 


single 


?■ 


The second occupancy effect comes from the initial 
rejection of background electrons, discussed in Section 


III C 2 where PMTs fired by background electron tracks 
are removed. If such an electron is close to a signal elec¬ 
tron in the RICH, the associated PMTs of the signal 
electron are also removed. The probability for this to 
happen is relatively small and increases with multiplic¬ 
ity. This loss is referred to as and it is estimated 

by monitoring the yield of e+e” pairs below 20 MeV/c^ 
before and after erasing the PMTs for each centrality 
bin. This mass region is dominated by Dalitz decays and 
7 conversions and provides a clean electron pair sample 
with a signal-to-background ratio of ^200 even for the 
most central events. Using these efficiency losses, 
can be expressed as: 


mult _ embed ^ TPMT 

^pair ^pair ’ ^pair 


( 22 ) 


Table 


IV 


for the five 


summarizes the values of and 

pail pail 

centrality bins. 



FIG. 22. (Color online) Pair efficiency correction for the pair 
Pt range between 0.8 and 1.0 GeV/c for each centrality bin. 
This represents the total efficiency including the elD selection 
cuts based on neural networks, losses in the acceptance due 
to detector inactive areas, losses induced by the pair cuts and 
occupancy effects in the central arm detectors. 


G. Systematic Uncertainties 

The main systematic uncertainties on the corrected 
data arise from uncertainties on the electron identifica¬ 
tion, the acceptance and the background subtraction. 
They are discussed in detail below and summarized in 
Table These uncertainties move all data points in the 
same direction but not by the same factor 


1. Systematic uncertainty on electron identification and 
occupancy effects 

As described in Section [III C| electron identification is 
achieved using three neural networks. Different thresh¬ 
old cuts for the neural networks result in different elec¬ 
tron identiheation efficiency and occupancy effects. The 
thresholds in the neural networks are varied by ±20% 
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TABLE V. Summary of systematic uncertainties assigned to the corrected data for MB collisions. 


Component 

Mass range 

Systematic uncertainty 


elD -I- occupancy effects 
Acceptance (time) 

Acceptance (MC) 

Combinatorial background 

0-5 GeV/c^ 

±4% 

±8% 

±4% 

±25% (mee = 0.6 GeV/c^) 


Residual yield 

0-0.08 GeV/c^ 

-5% (mee = 0.08 GeV/c^) 


Residual yield 

1-5 GeV/c^ 

-15% (mee = 1 GeV/c^) 



around the selected values and the variations of the elec¬ 
tron pair yield in the mass region rriee < 150 MeV/c^, 
after applying the efficiency correction, are used to as¬ 
sess the systematic uncertainty of electron identification 
and occupancy effects. 

By changing the thresholds by ±20% the raw electron 
pair yield changes by about ±50%. However, once the 
corresponding efficiency corrections are applied, the vari¬ 
ations are below 4% for all the centrality bins. Based on 
these results, we assign a ±4% systematic uncertainty on 
the electron identification. 


2. Systematic uncertainty on the acceptance 


We consider two sources of systematic uncertainties on 
the acceptance: variations of the pair acceptance vs time 
and variations of the pair acceptance between data and 
MC simulations. 

The pair acceptance systematic uncertainty vs time is 
studied by considering the variations of the number of 
electron pairs per event for each run group. The weighted 
average of the rms of the number of electrons per event 
in the five run groups is found to be 8% and it is taken 
as the systematic uncertainty of the acceptance variation 
over time. 

The systematic uncertainty on the data vs MC pair 
acceptance is studied by comparing the reconstructed 
yield in data and simulations. In data we select recon¬ 
structed pairs with uiee < 100 MeV jc ?, after subtracting 
the combinatorial and correlated components of the back¬ 
ground, using data from one of the run groups. In the 
MC simulations we use reconstructed pairs in the same 
mass range from 7r° Dalitz decays applying the fiducial 
cuts for the corresponding run group. The entire detec¬ 
tor is divided into four sectors. Data and MC simulations 
are normalized in one sector. The variations of the yield 
ratios between data and MC simulations in the other sec¬ 
tors ranges between 1% and 8%. The weighted average 
of these variations is found to be 4% and it is taken as 
the systematic uncertainty of the acceptance agreement 
between data and MC simulations. 


3. Systematic uncertainty on the background subtraction 


We consider two sources of systematic uncertainties on 
the background subtraction: 

(i) Uncertainty on the combinatorial background sub¬ 
traction. It is primarily due to the uncertainty in the 
normalization factor, and the latter is determined by 
the statistics in the normal ization window, namely by 
1/\/Nls (see Section IIIE6|). This translates into a rel¬ 


ative uncertainty of the signal SS/S = l/y/N^s x B/S. 
The ratio B/S depends both on mass and centrality. In 
Table [v| we quote the uncertainty at rriee = 0.6 GeV/c^ 
which represents the worst case in mass, for MB events. 
The centrality dependence results in variations of the or¬ 
der of 15% from the MB values. 

(ii) In the ideal case, the like-sign residual yield, i.e. 
the like-sign yield after subtracting all the background 
sources, should be zero. In practice it is not. As shown 
in Figs. [T6]and[l7l there is a small residual yield. In this 
analysis, we assume that any residual yield is entirely due 
to unsubtracted background, and we take it as an addi¬ 
tional source of systematic uncertainty, after transform¬ 
ing it into unlike-sign residual yield via the acceptance 
correction factor a. This uncertainty takes into account 
any possible discrepancy in shape or magnitude of the 
various subtracted sources of background. The factor a 
accounts for the different acceptance of the PHENIX de¬ 
tector for like and unlike sign pairs. It is calculated as a 
function of pair mass and pair pT using the mixed event 
background as: 


\.{m,pT) = 


MIX^ {m,pT) 


MIXj^+{m,pT) ± MIX _ (m,pT) 


(23) 


Figure 23 panels (a)-(e) show a times the like-sign 
residual yield divided by the sum of all unlike-sign back¬ 
ground sources as a function of mass for the five central¬ 
ity bins, which represent the relative residual background 
yield in the unlike sign mass spectrum. The mass regions 
ruee < 0.08 GeV/c^, 0.2 GeV/c^ < niee < 1.0 GeV/c^ 
and rriee > 1 GeV/c^ are fitted to a constant to quan¬ 
tify the magnitude of the residual unlike-sign yield. The 
fit results are also shown. Figure 23 panels (f)-(j) show 
zoomed views in the vertical axis for the 0.2-1 GeV/c^ 
mass range. The fits in the mass region rriee = 0.2- 
1.0 GeV/c^ give results that are consistent with zero for 
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FIG. 23. (Color online) (a-e) Unlike-sign residual background yield derived from the like-sign residual yield, obtained after 
subtracting all background sources, via the acceptance correction factor a (see text). The legend and the dashed lines show 
the results of constant fits below 80 MeV/c^ and above 1 GeV/c^. (f-j) Zoomed views in the vertical axis for the 0.2-1 GeV/c^ 
mass range. 


all centrality bins. For the other two mass ranges, the 
residual yields are considered as sources of systematic 
uncertainties if their signihcance is larger than 2cr. 


The total systematic uncertainty in the background 
subtraction is obtained as the quadratic sum of the 
systematic uncertainties due to the combinatorial back¬ 
ground subtraction and the residual yield. Both contri¬ 
butions are listed in Table IVl for MB collisions. It is worth 
noting that the systematic uncertainty of the background 
subtraction is much lower than the required accuracy to 
measure a signal with the S/B values shown in Section 
IIIIE71 


H. Cross checks 


A second independent analysis was performed as a 
cross check. The key features of the second analysis are 
discussed here. A more detailed description is given in 
Appendix The second analysis is similar to the anal¬ 
ysis described in Ref. [23], but it makes use of the HBD 
and includes all the important improvements developed 
in this work. In particular, it makes use of the time-of- 
flight information for better hadron rejection, implements 
the shape distortion of the mix ed event background due 
to elliptic flow (Section III E 2), subtracts the correlated 
electron-hadron background (Section III E 5), and explic¬ 
itly considers the away-side jet-pair component in the 
background subtraction (Section [ill E4). 


Important elements of the independent analysis are dif- 
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ferent from those of the main analysis. The most signifi¬ 
cant differences are: (i) The HBD underlying event sub¬ 
traction is done using the average charge in the vicinity of 
a track as opposed to the average charge in a module as 
used in the main analysis, (ii) Electron identification is 
achieved by a sequence of independent one-dimensional 
cuts on each of the electron identification variables in¬ 
stead of the neural network approach, (iii) The normal¬ 
ization of each background source is determined from a 
fit to the like-sign spectra, in contrast to the main analy¬ 
sis where all the correlated background sources are abso¬ 
lutely normalized and only the combinatorial background 
is normalized to the like sign spectra. 

The second analysis results in a factor of two smaller 
signal-to-background ratio and a 10 % reduction in purity 
of the electron sample in central collisions. However, once 
corrected for efficiency, the results of the second analysis 
are consistent within uncertainties with those obtained 
with the main analysis described in this section. 


IV. COCKTAIL OF HADRONIC SOURCES 

In this section we describe the procedures used to cal¬ 
culate the expected dielectron yield from hadronic de¬ 
cays, commonly referred to as the hadronic cocktail, that 
will be compared to the experimental results in Section 
The known e“''e“ sources are calculated using the EX¬ 
ODUS, PYTHIA and MC@NLO event generators. EXODUS is 
a phenomenological generator that simulates phase space 
distributions of the relevant electron sources and their de¬ 
cays [SO]. It generates the photonic sources, i.e. Dalitz 
decays of light neutral mesons: 7r°, 77, 77 ' —>■ and 

ui —>■ e^e~TT^ and the nonphotonic sources, i.e. dielectron 
decays of mesons: p, w, <j), J/t/i —>• e+e“. PYTHIA [55] 
and MC@NLO [STJ1^ are used to generate the correlated 
pairs from semi-leptonic decays of heavy flavor (charm 
and bottom) mesons. The hadrons are assumed to have 
uniform pseudorapidity density within |? 7 | <0.35 and uni¬ 
form azimuthal distribution in 27r. Once generated, the 
sources are filtered through the ideal acceptance of the 
PHENIX detector and smeared with the detector res¬ 
olution for comparison to the measured invariant mass 
spectrum. 


A. Neutral pions 


The dominant electron source as well as the funda¬ 
mental input for EXODUS is 7r°. The shape of the 7r° px 
distribution is parameterized as: 


^ _ 1 __ 

(Pp {e-^P-r-bpl 


(24) 


The parameters, a, 6 , po and n, are obtained by a simulta¬ 
neous fit of the PHENIX published results for 7 r° 
and charged pions [55j . The resulting fit parameters are 


shown in Table [Vl] for the five centrality bins of this anal¬ 
ysis. The absolute magnitude of the 7 r° rapidity density, 
dN^o/dy, is obtai ned by fitting the cocktail to the data 
(see Section IVD). 


TABLE VI. Fit parameters derived from the tt® and charged 
pion pt distributions |53H55 | for different centralities using 
Eq. (|2t|. 


Parameter 

0%-10% 

10%-20% 

0 

6 ? 

0 

CN 

40%-60% 

60%-92% 

a [(GeV/c)-b 

0.57 

0.53 

0.43 

0.36 

0.33 

b [(GeV/c)- = ] 

0.19 

0.16 

0.11 

0.13 

0.088 

po [GeV/c] 

0.74 

0.75 

0.79 

0.76 

0.74 

n 

8.4 

8.3 

8.5 

8.4 

8.4 


B. Other mesons 


The Pt distributions of other light mesons are based 
on the parametrization of the pion spectrum assuming 
rriT scaling [53], i.e. Eq. (241 is used withpr replaced by 

— mio ■ This assumption reproduces well 




Pt 


m. 


meson 


the measured light meson px distributions in Au-|-Au col¬ 
lisions as demonstrated in |23j . The absolute normaliza¬ 
tion for each meson is provided by the ratio of the meson 
to 7 r° invariant yield at high px {px > 5 GeV/c ). W e use 
the values from Ref. [55j, summarized in Table [vTT 


TABLE VII. Meson to tt® ratio at high px {px > 5 GeV/c) 
obtained from PHENIX data in p+p collisions [44] . 


77 / 77 ° 

p/7r° 

57 / 77 ° 

77777 ° 

<//77° 

0.48 

1.0 

0.90 

0.25 

0.40 


The values were obtained from p+p collisions and are 
taken to be valid for Au-|-Au collisions because at high px 
the suppression of all mesons is found to be very similar 
to the TT® suppression and consequently the meson/ 7 r*^ ra¬ 
tios in Au-f Au collisions remain unchanged with respect 
to the ratios in p+p collisions [56H58] . 

For the px distribution of the J/tp we use the neu¬ 
tral pion Px spectrum measured in p+p collisions m, 
assuming mx scaling. Detector effects on the J/ip line 
shape are taken into account by passing the decay e“''e“ 
through a GEANT simulation of the PHENIX detector. 
The resulting pr integrated invariant e“'"e“ mass distri¬ 
bution is then normalized to the measured cross section 
in p+p collisions [23] and scaled to Au-|-Au collisions by 
the corresponding (fVcoii) and the measured Raa for each 
centrality bin [35] . 


C. Open heavy flavor 

The correlated e“'"e“ yield from open heavy flavor de¬ 
cays is simulated using two different p+p event gener- 
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ators, PYTHIA and mc@NLO, and measured cc and bb 
production cross sections. 

PYTHIA simulations are used to calculate gluon fu¬ 
sion, the dominant process for heavy-quark production, 
in leading-order perturbative QCD. Specifically, we use 
PYTHiA-6 [50] and CTEQ 5 L as input parton distribution 
functions. The MC@nlo package (vers. 4.03) jOU [50] 
is a next-to-leading order simulation that generates hard 
scattering events. These events are subsequently fed to 
HERWIG (vers. 6.520) jOT] for fragmentation in vacuum. 

We use the cc- and 66-production cross sections mea¬ 
sured by PHENIX [00], by fitting the event generator 
(pYTHiA or MC@NLO) Output to the measured dielec¬ 
tron mass spectrum in d-|-Au collisions for me+e- > 
1.15 GeV/c^. These cross sections were scaled by the 
average number of d-|-Au binary collisions ((A^coii)) to 
give the p+p equivalent cross section. For bb, both gen¬ 
erators gave within uncertainties the same result for the 
cross section extrapolated to zero invariant mass |62] : 


dy 


y=o 


1.36 ± 0.32(stat) ± 0.44(syst) ^b 


(25) 


The cc cross section strongly depends on the event 
generator. The mc@NLO yields the following cross sec¬ 
tion [52]: 


daP? 


dy 


v=o 


= 287 ± 29(stat) ± lOO(syst) ^b (26) 


whereas PYTHIA gives: 
daP? 


dy 


y=o 


= 106 ± 9(stat) ± 33(syst) ^b (27) 


This cross section, derived from e+e“ data in d-|-Au col¬ 
lisions, is consistent within uncertainties with the cross 
section derived from measurements of single electrons 
from semileptonic decays of heavy flavor mesons in p+p 
collisions, extrapolated to pt = 0 GeV/c using pythia 
simulations jH]. mc@NLO was not used to derive the 
heavy flavor cross section from measurements of single 
electrons. 


The two results, Eqs. (26) and (27), although con¬ 
sistent within ^1.2 a, yield central values which dif¬ 
fer by a factor of ^2.5. This difference comes mainly 
from the extrapolation of the dilepton yield from ruee 
> 1.15 GeV/c^ to TOee = 0 GeV/c^, as illustrated in 
Fig. Figure [Mj also shows an absolute comparison of 
the PYTHIA and MC@NLO dielectron invariant yields from 
correlated heavy flavor meson decays in MB Au-|-Au colli¬ 
sions, obtained by IVcoii scaling of the p+p cross sections 
quoted in Eqs. (26) and (27). At high masses, mge > 


2 We use PYTHlA-6 Eq] with the following parameters MSEL[cc]=4 
or MSEL[66]=5, MSTP(91)=1, PARP(91)=1.5, MSTP(33)=1, 
PARP(31)=1.0, MSTP{32)=4, PMAS(4)=1.25, PMAS(5)=4.1. 



FIG. 24. (Color online) Comparison of the invariant dielec¬ 
tron yield from correlated heavy flavor meson decays for MB 
Au-|-Au collisions calculated with pythia (solid line) and 
MC@NLO (dashed line) using the /dy cross sections of 
106 ph and 287 ph, respectively [62] , scaled by (Acoii). 


1.15 GeV/c^, both generators give by construction the 
same yield, with a very small difference in shape. How¬ 
ever, at low masses there is a large discrepancy in the 
absolute yield. 

The d-l-Au (as well as the p+p) inclusive dilepton yield 
is not very sensitive to this variation of the cross section 
because the large effect at low masses is diluted by the 
contributions from light meson decays. The situation is 
quite different in Au-f Au collisions. The yield from light 
meson decays scales approximately with A^part, whereas 
the contribution from heavy flavor scales with iVcoii mak¬ 
ing the latter dominant at low-masses in central colli¬ 
sions. The choice of the generator used to simulate the cc 
contribution will therefore affect the total cocktail yield 
at low masses and will influence the interpretation of the 
Au-|-Au data in terms of an excess with respect to the 
cocktail. The results will be presented in the next sec¬ 
tion using PYTHIA for an easier comparison with previ¬ 
ously published results but both generators, pythia and 
MC< 0 lNLO, will be considered in the discussion. 


D. Cocktail normalization 


In the present analysis we use the precisely measured 
e“'"e“ data at low masses to derive the normalization of 
the cocktail of hadronic sources. In the restricted phase 
space defined by rriee < 0.1 GeV/c^ and prlmee > 5 
the inclusive e“'"e“ yield is dominated by tt^ Dalitz de¬ 
cays with a small contribution of direct virtual photons 
and an even smaller contribution of rj Dalitz decays. To 
a very good approximation the mass spectrum of these 
three sources has a l/ruee dependence and their relative 
magnitude is well known. The ratio of direct photons to 
7 r° is known from PHENIX measurements [531 [M] and the 
ratio of r] to 7r° can be easily obtained from the PHENIX 























26 


measurement at high pT |58j and the rriT scaling as de¬ 
scribed in Section |IVB By fitting the cocktail-l-direct 
virtual photons to the data in the restricted phase space 
defined above, one obtains the rapidity density dN^^o/dy 
that determines the normalization of the cocktail. The 
values are found to be consistent with measurements of 
neutral and charged pions |53H55j within the systematic 
uncertainties of cocktail and data. 

Alternatively, the cocktail can be absolutely normal¬ 
ized using the 7r° rapidity density dN^j^o /dy derived from 
these measurements as done in Ref. [23]. The cocktails 
obtained with these two procedures are compared in Fig. 
2^ The results differ at masses niee < 100 MeV/c^ by 


about 25% which is approximately the contribution of 
the virtual direct photons. However, for the mass range 
of interest, that is typically 0.3-0.76 GeV/c^, the differ¬ 
ence is smaller and amounts to only 15%. In this mass 
range, the yield is dominated by the contributions from 
correlated heavy flavor decays and changing dN^j^o /dy by 
^25% has a minor effect on the inclusive e+e“ yield. At 
even higher masses, niee > 1 GeV/c^, the two procedures 
yield exactly the same results. The present procedure is 
adopted to be consistent with the known contribution of 
internal conversion. 



FIG. 25. (Color online) Cocktail of hadronic sources for the 
2010 run with normalization provided by fitting to the present 
e^e“ invariant yield at masses niee < 0.1 GeV/c^ (black line) 
or with absolute normalization to the 7r° rapidity density de¬ 
rived from measurements of neutral and charged pions |53H55| 
(dashed line). 


(ii) Direct photon: The systematic uncertainties in the 
direct photon dN/dy are taken from Ref. [64]. They 
range from ±24% to ±70% from central to peripheral 
collisions, respectively. 


(hi) Open heavy flavor (cc, bb): We use the systematic 
uncertai nties of t he o pen heavy fl avor cross sections given 
in Eqs. ( [2^ or (27) for cc and (25) for bb, taken from 
Ref. |52|. The (Ncoii) systematic uncertainties shown in 
Table |ll| are added in quadrature when the p+p cross 
sections are scaled to Au±Au collisions. 


(iv) J/i/j: The systematic uncertainty of the J/V' cross 
section in p+p collisions is estimated to be ±14% [65]. 
The systematic uncertainties in (A^coii) and J/ip Raa are 
added in quadrature. The Raa uncertainties are taken 
from Ref. [SH], ranging from ±22% to ±35% depending 
on centrality. 


A summary of the cocktail systematic uncertainties is 
presented graphically in Fig. 26 which shows the system¬ 
atic uncertainty of each cocktail component together with 
the total cocktail systematic uncertainty, determined as 
their quadratic sum. 



FIG. 26. (Color online) Systematic uncertainties assigned 
to each cocktail component and the total cocktail systematic 
uncertainty for MB events. 


E. Systematic uncertainties on the cocktail 

The systematic uncertainties of the cocktail ingredi¬ 
ents are estimated and propagated to determine the total 
cocktail systematic uncertainty. The following uncertain¬ 
ties are considered: 

(i) Light meson to 7r° ratio: We adopt the same sys¬ 
tematic uncertainties used in Ref. [53] , namely ±30% for 
rj, uj and (j), ±33% for p and ±100% for y'. 


F. The Au+Au Cocktail 

The cocktail, calculated as described above, using the 
PYTHIA generator for the open heavy flavor contribu¬ 
tions, is presented in Fig. [^ for MB Au±Au collisions 
together with the individual components of the cock¬ 
tail. For comparison. Fig. [^ also shows the total cock¬ 
tail using MC@NLO for the open heavy flavor contribu¬ 
tions. The differences discussed above in Section IIV Cl 
are clearly reflected in this comparison. 
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FIG. 27. (Color online) Cocktail of hadronic sources for the 
2010 run (black solid line) using the pythia generator for the 
open heavy flavor contributions. The individual components 
of the cocktail are also shown. For comparison, the total 
cocktail using mc@nlo is shown (black dashed line). 


V. RESULTS AND DISCUSSION 
A. Invariant mass spectra 



FIG. 28. (Color online) Invariant mass spectrum of e^e“ 
pairs in MB Au+Au collisions within the PHENIX acceptance 
compared to the cocktail of expected decays. 


0.76 GeV/c^, below the p meson mass, that is the most 
sensitive region to in-medium effects. 


Figure 28 shows the invariant mass spectrum of e“'"e“ 
pairs within the PHENIX acceptance (as defined in Sec¬ 
tion IIEl) for MB Au-|-Au collisions. The spectra are 


subject to a pt cut of 0.2 GeV/c on the single electron 
tracks and to a 100 mrad cut on the pair opening an¬ 
gle. Statistical and systematic uncertainties on the data 
points are shown separately by vertical bars and boxes, 
respectively. Figure also compares the measured spec¬ 
trum to the cocktail of expected e+e“ sources, where 
PYTHIA is used to calculate the correlated pairs from 
heavy flavor decays. The individual contributions to the 
cocktail are shown in the figure. 

See Section |IV] for details about the cocktail calcula¬ 
tion. The total systematic uncertainty of the cocktail is 
shown by the yellow band. The bottom panel shows the 
ratio of data to cocktail. 

Figure 29 shows the invariant mass spectra of e“'"e“ 


pairs for the five centrality bins analyzed in this work, 
compared to the cocktail. 

For a more detailed discussion of the centrality and 
transverse momentum dependencies of the dielectron 
yield, we consider three mass regions: 

(a) the mass region ruee < 0.10 GeV/c^ that is 
dominated by the 7r° Dalitz decay. 


(b) the low-mass region (LMR), 0.30 < rUe 


< 


(c) the intermediate-mass region (IMR), 1.2 < mee < 
2.8 GeV/c^, that is dominated by the correlated pairs 
from the semi-leptonic decays of charm and bottom 
mesons. 

Figurej^shows the pair pT distribution for these three 
mass intervals in MB collisions. In the following sections 
we discuss the results in these three mass intervals. 


B. 7r° Dalitz region 


The mass region Wee < 0.10 GeV/c^ is dominated by 
the TT*^ Dalitz decay with a small contribution of direct 
virtual photons of ^20% and an even smaller contribu¬ 
tion of the rj Dalitz decay of ~10%. We discuss here only 
the shape of the pt distribution because the integrated 
dielectron yield in this mass interval was used to normal¬ 
ize the cocktail for the five centrality bins as described in 
Section |IV] Figure [30| compares the measured dielectron 
Pt distribution for MB collisions in this mass interval to 
the Pt distribution of the hadronic cocktail that uses the 
parametrization for the 7r° and rj mesons [Eq. (24|]. The 
agreement between the two distributions, in shape and 
magnitude, is very good when adding the measured yield 
of direct virtual photons. 
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in Au+Au collisions within the PHENIX acceptance for the 
various centrality bins. The lines represent the total expected 
yield from all the sources indicated in Fig. |28| 


C. Low-mass region (LMR) 


In the LMR, the yield is expected to be saturated by 
the lig ht m esons ( 77 , p and w) and the cc contribution. 
Figure 28 shows an enhancement of e“'"e“ pairs with re¬ 
spect to the cocktail in MB collisions. The enhancement 
develops with centrality as shown in Fig. and it ap¬ 
pears to be distributed over the whole pr range covered 
by the measurement, as can be seen in Fig. We 

quantify the effect by the enhancement factor defined as 
the ratio of the measured over expected dilepton yield 
integrated in the LMR. As discussed in Section [TV C| the 
cocktail yield in this mass region depends on the gen¬ 
erator, PYTHIA or MC@NLO, used to calculate the open 
heavy flavor contribution. The enhancement factors ob¬ 
tained with PYTHIA are shown as a function of centrality 
in Fig. I^and they are listed in in Table VIII for the two 
cases. The enhancement factors are approximately 40% 
higher when pythia is used to calculate the open heavy 
flavor contribution instead of mc@NLO. 



FIG. 30. (Color online) MB invariant pr distributions for 
three mass windows as indicated in the legend. The solid 
lines represent the expected pt distributions of the hadronic 
cocktail and the shadowed bands around the lines represent 
the cocktail systematic uncertainties. The dotted lines include 
the contribution from direct photons in the phase space region 
where they can reliably be calculated, i.e. pr/mee > 5. 


TABLE VIII. Enhancement factors, defined as the ratio of 
measured over expected dilepton yield in the mass region rriee 
= 0.30-0.76 GeV/c^, for the five centrality bins and for MB. 
The enhancement factors are quoted separately for the two 
cases where the correlated yield from cc decays is calculated 
with pythia or mc@nlo. The imodel uncertainties repre¬ 
sent the cocktail systematic uncertainties. 


Gentrality 

Enhancement factor ±stat ±syst Tmodel 


MC@NLO cc 

PYTHIA CC 

MB 

1.7 ±0.3 ±0.3 ±0.2 

2.3 ± 0.4 ± 0.4 ± 0.2 

0%-10% 

2.3 ±0.7 ±0.5 ±0.2 

3.2 ± 1.0 ±0.7±0.2 

10%-20% 

1.3 ±0.4 ±0.5 ±0.2 

1.8 ±0.6 ±0.7 ±0.2 

20%-40% 

1.4 ±0.2 ±0.3 ±0.2 

1.8 ±0.3 ±0.4 ±0.2 

40%-60% 

1.2 ±0.2 ±0.3 ±0.2 

1.6 ±0.2 ±0.4 ±0.2 

60%-92% 

1.0 ±0.1 ±0.2 ±0.2 

1.4 ±0.2 ±0.3 ±0.2 
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FIG. 31. (Color online) Data to cocktail (using pythia for 
heavy flavor contribution) ratio in the LMR versus centrality. 
The shaded band around one represents the cocktail system¬ 
atic uncertainty. 


• Pair cuts: Loose pair cuts were applied in Ref. [23] 
compared to the cuts used in this analysis. The 
cuts used in Ref. [23] are found to leave a sizable 
amount of detector induced correlation in the mass 
region ruee = 0.4-0.6 GeV/c^. 

• Flow: As demonstrated in Section IIIIE2I the col¬ 
lective flow that is inherent to nuclear collisions, 
affects the shape of the combinatorial component 
of the background and violates the square root re¬ 
lation [Eq. ®]- These two effects were not taken 
into account in the data analysis of Ref. [231 ■ 

• Electron-hadron pairs: As shown in Section pII E 5[ 
the e-h pairs originate in the central arm detectors 
and in particular in the RICH detector. This source 
of correlated pairs was not considered in [23] . 


1. Comparison to previous PHENIX results 

The enhancement factors qnoted above are sig¬ 
nificantly smaller than those previously reported by 
PHENIX [23] in the same Au-|-Au collision system at the 
same energy of -./s^ = 200 GeV. There are a number of 
significant differences, both qualitative and quantitative, 
between the two analyses: 

• Hadron contamination: The purity of the electron 
sample is very different in the two cases. In [23] the 
hadron contamination was 30% in central Au-|-Au 
collisions, whereas in the present analysis, the HBD 
enabled this contamination to be reduced to less 
than 5% at all centralities. 

• Signal sensitivity: The signal sensitivity is usu¬ 
ally quantified by the signal to background S/B 
ratio. The S/B values displayed in Eig. [I^ are 
similar to those quoted in Ref. [23] • This is how¬ 
ever, a misleading comparison, because in a situ¬ 
ation of subpercent S/B ratio, the magnitude of 
S critically depends on the accuracy of the back¬ 
ground subtraction. A better way to assess the 
sensitivity of the measurement is provided by the 
cocktail/background, C/B, ratio. From the sig¬ 
nal/background ratio and the enhancement factors 
quoted in Ref. [23] j we estimate an average value of 
C/B over the mass range Wee = 0.15-0.75 GeV/c^ 
of ~l/600 in MB collisions. In the present analysis 
the same ratio is found to be ~l/250. In addition 
to that, one should take into account that in the 
2010 run with the H— field configuration there is 
a larger track acceptance of ^20%. This rough es¬ 
timate indicates that at the same multiplicity the 
signal sensitivity in the present analysis is larger by 
a factor of ^3.5 compared to the previous one. 


• Away-side jet component: The away-side jet com¬ 
ponent of the correlated background was found to 
be negligible in |23j and only the near-side jet com¬ 
ponent was considered. In the present analysis, 
both components are absolutely calculated. The 
away-side component is indeed relatively small but 
both components are considered and subtracted. 

• Background subtraction procedure: In Ref. [23] . 
the shapes of the three components of the back¬ 
ground (combinatorial background, cross pairs and 
near-side jet) were calculated whereas their abso¬ 
lute scales were obtained by fitting to the like-sign 
spectra. In the present analysis, all components 
of the correlated background (cross pairs, jet pairs 
and electron-hadron pairs) are calculated and sub¬ 
tracted in absolute terms. There is only one free pa¬ 
rameter in the background subtraction procedure, 
namely the normalization factor of the combinato¬ 
rial background. 

In conclusion, we do not confirm our previous report of 
a large excess seen in the LMR [23] . The differences listed 
above affect the yield in the mass region where the excess 
was reported but not always in the same direction. For 
example, the loose pair cuts lead to under subtraction of 
the background whereas neglecting the flow modulation, 
has the opposite effect namely it leads to over subtraction 
in the mass region where the excess was observed. These 
differences also do not affect the unlike-sign yield by a 
similar magnitude. The hadron contamination, the loose 
pair cuts and the electron-hadron pairs are the most sig¬ 
nificant ones in this respect. Taking all the differences to¬ 
gether, the present analysis is much improved compared 
to the previous one and we thus consider the previous 
result on the low-mass excess to be superseded by the 
results presented here. 
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2. Comparison to STAR results 

Recently, STAR published results on e+e“ production 
in Au+Au collisions at = 200 GeV [Ml |H7] • In 

the same mass range of Wee = 0.30-0.76 GeV/c^, STAR 
observes an excess of dielectrons and quotes a value of 
1.77±0.11®*“‘ ± 0.24®2^®* ± 0.33'"°'^®' in MB collisions, for 
the ratio of the dielectron yield to the hadronic cock¬ 
tail excluding the p meson contribution. There are two 
factors that should be taken into account when compar¬ 
ing the STAR results with those quoted in Table |VIII| 
First, excluding the p contribution results in an increase 
of about 10% of the data to cocktail ratio. Second, STAR 
uses PYTHIA with a charm cross section dacc/dy = 171 
±26 ph |M] which is between the PHENIX cross sections 
quoted in Section |TV] for pythia and mc@NLO. Taking 
those two differences into account, as well as the experi¬ 
mental uncertainties, we find that the results of the two 
experiments are consistent in the LMR. The centrality 
and pt dependencies of the enhancement reported in m 
are also consistent with our results. 
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FIG. 32. (Color online) Data to cocktail ratio in the IMR 
versus centrality. The cocktail uses pythia for the cc con¬ 
tribution (left scale) or random cc contribution (right scale). 
The shaded band represents the pythia cocktail systematic 
uncertainty. The same uncertainty applies also to the random 
cc cocktail. 


TABLE IX. Enhancement factors, defined as the ratio of mea¬ 
sured to expected dilepton yield in the mass region rriee = 
1.2-2. 8 GeV/c^, calculated using pythia for the five central¬ 
ity bins and for minimum bias. The last line gives the en¬ 
hancement factor assuming random correlation (see text). 



Gentrality 

Enh. factor ±stat ±syst ±model 

PYTHIA CC 

0 %-10% 

1.3 ±0.7 ±0.2 ±0.3 


10 %-20% 

1.8 ±0.5 ±0.3 ±0.3 


20%-40% 

1.8 ±0.2 ±0.3 


40%-60% 

1.1 ±0.2 ±0.1 ±0.3 


60%-92% 

1.0 ±0.2 ±0.1 ±0.3 


MB 

1.5 ±0.3 ±0.2 ±0.3 

MB (random 

cc) 2.5 ±0.5 ±0.3 ±0.3 



D. Intermediate-mass region (IMR) 

The IMR is dominated by correlated pairs from the 
semi-leptonic decays of DD mesons, with a small contri¬ 
bution from BB mesons and an even smaller contribution 
from Drell Yan. The latter is neglected in the cocktail 
calculation. This mass interval is singled out by theory 
as the most sensitive window to identify the thermal ra¬ 
diation of the QGP in the dilepton spectrum [551169j . 

The results displayed in Figs. [28] and [29| show a small 
enhancement of dileptons in the IMR with respect to the 
yield from cc decays calculated using pythia. The en¬ 
hancement factors are shown in Fig. [^ as a function 
of centrality and the values are listed in Table |IX| The 
results are consistent with those of Ref. [53] within the 
large experimental uncertainties of the latter. There is 
very little difference in the dilepton yield in this mass 
interval if MC@NLO is used instead of pythia, as demon¬ 
strated in Fig. [^ The shapes are very similar and the 
integral yields in the IMR differ by less than 10% in the 
two cases. 

Using PYTHIA, the enhancement factor in MB events 
is ~1 standard deviation away from unity. However, the 
data to cocktail comparison discussed above, represents 
an extreme case in which it is assumed that the cor¬ 
relations between the cc pairs in Au±Au collisions are 
the same as in p+p collisions. It is however, well known 
that heavy flavor quarks exhibit energy loss and collec¬ 
tive flow in the medium formed in Au±Au collisions, as 
manifested for example in measurements of single elec¬ 
trons [331170] . This should affect the correlation between 
the e+e“ pairs from cc decays. Lacking a suitable gen¬ 
erator to model this effect, we consider also the opposite 
extreme approach in which we assume that the pair is 
totally decorrelated. The invariant mass is calculated us¬ 


ing two electrons randomly selected from the measured 
Pt distribution of single electrons from heavy flavor de¬ 
cays [33], with uniform distributions in pseudorapidity 
and azimuthal angle. The pair is Hltered through the 
ideal PHENIX acceptance and the integral is normalized 
to the calculated pythia yield from cc decays. This ex¬ 
treme case results in a softer mass distribution in the 
IMR as can be seen in Fig. [^ 

There is a small yield depletion at high masses com¬ 
pensated by a higher yield at low masses. The integral in 
the IMR is lower resulting in enhancement factors that 
are ~70% larger compared to those derived from pythia. 
The enhancement factor in MB collisions is quoted in the 
last line of Table [TX] and the centrality dependence is seen 
by comparing the data points to the dot-dashed line in 
Fig.[3g 
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FIG. 33. (Color online) Invariant mass spectrum of e^e“ 
pairs in MB Au+Au collisions within the PHENIX acceptance 
compared to the cocktail of expected decays when the cc decay 
component is calculated assuming no correlation between the 
c and c. 


E. Comparison to theory 

In this section we compare our results to the model 
originally developed by Rapp and Wambach [7I1II2]. The 
model uses an effective Lagrangian and a many body ap¬ 
proach to compute the electromagnetic spectral function 
which is the main factor in the calculation of the dilepton 
production rates. In the LMR, the spectral function is 
saturated via vector meson dominance, by the light vec¬ 
tor mesons, in particular the p meson, whereas at larger 
masses it is dominated by multipion states or equiva¬ 
lently, via quark-hadron duality, by qq annihilation. The 
dilepton yields are obtained by an appropriate integra¬ 
tion of the thermal rates over the space-time evolution 
of the fireball. This model was very successful in repro¬ 
ducing the low-mass dilepton enhancement discovered at 
SPS by the CERES experiment and later further studied 
by the NA60 experiment. In the comparison below, we 
use an improved version of the model that incorporates 
recent developments, a nonperturbative QGP equation 
of state and QGP emission rates, i.e. qq annihilation at 
temperatures higher than the critical temperature, both 
based on lattice QCD [73]. It is important to note that 
this updated version preserves the agreement with the 
SPS data and also reproduces the RHIC results from 
STAR. 

Figures and compare the invariant mass spec¬ 
trum and the LMR pair px distribution with the model 


calculations for MB collisions m- The main compo¬ 
nents, in-medium p broadening, QGP thermal radiation 
and cocktail excluding the p, together with their sum, 
are shown separately. 



FIG. 34. (Golor online) MB invariant mass spectrum com¬ 
pared to the model calculations of Rapp (solid line) | 74| . 
The main contributions, the in-medium p broadening (dot¬ 
ted line), the QGP thermal radiation (dot-dashed line) and 
the cocktail excluding the p (dashed line) are also shown. 



FIG. 35. (Golor online) Dielectron pt distribution in the LMR 
compared to model calculations (solid line) [73]. The main 
contributions, the in-medium p broadening (dotted line), the 
QGP thermal radiation (dot-dashed line) and the cocktail ex¬ 
cluding the p (dashed line) are also shown. 

In both figures the data are consistent with the cal¬ 
culations. Within this model, the enhancement in the 
LMR originates from the in-medium p broadening, i.e. 
the thermal radiation of the hadronic phase, with a very 
small contribution from the QGP. 

In the model, the centrality dependence of the ther¬ 
mal radiation is reasonably well described, within an 
uncertainty of ~I0%, by a power-law scaling of the 
charged particle rapidity density {dNch/dy)°‘, with a ~ 
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FIG. 36. (Color online) Centrality dependence of the dielec¬ 
tron excess, defined as (data — cocktail excluding p) compared 
to the thermal radiation from the hadronic (p broadening) and 
QGP phases from model calculations (dashed line) [74j . 


1.45 [73], very similar to the scaling of the thermal pho¬ 
ton yield |64l[69]. Within uncertainties, the present data 
are consistent with this scaling as illustrated in Fig. 
which also shows the centrality dependence of the excess, 
i.e. the data after subtracting the cocktail without the 
vacuum p, together with the expected power-law scaling 
(dashed line). 


VI. SUMMARY AND CONCLUSIONS 

PHENIX has measured invariant mass spectra, px dis¬ 
tributions and the centrality dependence of the e“'"e“ pair 
production in Au-I-Au collisions at = 200 GeV. 

The use of the HBD provided additional electron identifi¬ 
cation to the central arm detectors, additional hadron re¬ 
jection and increased rejection of the combinatorial back¬ 
ground. 

A new analysis procedure based on neural networks 
has been developed that combines in an efficient way the 
information from the HBD and the central arm detec¬ 
tors, RICH, TOF and EMCal. This results in three in¬ 
dependent parameters for electron identification, hadron 
rejection and close pair rejection, instead of the fourteen 
parameters of the four detectors involved in these tasks. 
A quantitative understanding of the total background at 
the subpercent level is achieved in the most central col¬ 
lisions. This is realized by a precise evaluation of all the 
background sources. The combinatorial background is 
determined by the event mixing technique together with 
an exact weighting procedure to take into account the 
flow effects that are inherent in the foreground events 
and cannot be reproduced in the mixed events. All the 
correlated background sources are calculated in absolute 
terms using simulations and published results. 

The results are compared with a cocktail of the known 
e“'"e“ sources. The contributions from light hadron de¬ 


cays that dominate the e“'"e“ yield at low masses rriee < 
1 GeV/c^, are determined using PHENIX measurements 
for pions and mx scaling for other mesons. The contri¬ 
butions from semileptonic decays of heavy flavor (charm 
and bottom) mesons are calculated with the pythia or 
MC@NLO generators using (Ncou) scaled p+p cross sec¬ 
tions. Both generators give very similar yields in the 
IMR. However, they predict very dissimilar results that 
differ from each other by a factor of ~2 in the LMR. 
Precise measurements of the charm cross section over the 
entire phase space are needed to resolve this discrepancy. 

A small enhancement of e“'"e“ is observed in the 
LMR with respect to the cocktail. The enhancement 
is distributed over the entire px range measured {px < 
5 GeV/c). It increases with centrality and amounts to 
2.3±0.4(stat)±0.4(syst)±0.2™°‘^®' for MB collisions when 
PYTHIA is used to calculate the open heavy flavor contri¬ 
bution. If instead mc@NLO is used, the enhancement fac¬ 
tors are ^40% smaller and for MB collisions it is found to 
be 1.7±0.3(stat)±0.3(syst)±0.2™°‘^®b The large enhance¬ 
ment of e“'"e“ pairs in the LMR previously reported by 
PHENIX, in Au-|-Au collisions at -^/s^ = 200 GeV [53| , 
is not confirmed by the results of the present improved 
analysis. In particular, the concentration of the excess at 
lowpT {pt < 1 GeV/c) is not observed here. The present 
results are consistent with those recently published by the 
STAR Collaboration [^ within the uncertainties of the 
two experiments. 

In the IMR, the results are compared with calcula¬ 
tions of the expected yield from the semileptonic decays 
of heavy flavor mesons in two extreme scenarios. In the 
first scenario, the heavy flavor contribution is calculated 
assuming that the correlations between the cc are the 
same in Au-I-Au as in p+p collisions, ignoring decorrela¬ 
tion effects produced by the interactions of heavy flavor 
quarks with the medium. A small enhancement is ob¬ 
served with respect to the yield predicted by pythia. It 
amounts to 1.5 ± 0.3(stat) ± 0.2(syst) ± 0.3“°'^®* for MB 
collisions. In the other scenario, the opposite extreme 
approach is adopted where the pair is assumed to be to¬ 
tally decorrelated. In this case, the enhancement factor 
becomes 2.5 ± 0.5(stat) ± 0.3(syst) ± 0.3“°'^®*. The re¬ 
ality is somewhere between these two extreme cases and 
we conclude that there is room in the data for a sig¬ 
nificant additional contribution, for example of thermal 
radiation, in the IMR. The nature of the IMR pairs will 
be studied with high statistics Au-I-Au data in 2014 data 
taking with the silicon vertex tracker (VTX) installed in 
PHENIX. 

The results in the LMR are compared to calculations 
based on the model originally developed by Rapp and 
Wambach Iniiig with subsequent improvements that in¬ 
corporate recent developments m- The model includes 
thermal radiation emission from the QGP phase {qq an¬ 
nihilation) as well as from the hadronic phase (mainly 
from the p meson copiously produced by pion annihi¬ 
lation, 7r''"7r“ —I p —I e^e“). The invariant mass and 
Px distributions as well as the centrality dependence are 
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well reproduced by the calculations. The enhancement 
observed in the LMR from SPS up to RHIC energies is 
thus consistently reproduced by a single model. Within 
this model, the enhancement originates from the melting 
of the p meson resonance as the system approaches chiral 
symmetry restoration. 
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Appendix A: Introducing Flow in the Mixed Events 


Pfg{^(P) 
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(A2) 

dpidp2P(Pi - df)P(p2 - 


/ 7r/2 ^TT 

d'k / dpiP(pi — df)P(pi + Ap — df) 
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Replacing P{p — d>) hy its expression in ( |Al[ ) allows 
one to write Pfg as the sum of four integrals: 


Pfg{AP) = - f ^ ddf f dPi{A + B + C + DOAS) 

J —7r/2 J —TT 


A = e(pi)e(pi + Ap) (A4) 

B — 2v2e{pi)e{pi + Ap) cos 2{pi — Jt) (AS) 
C = 2v2e{pi)^{pi + Ap) cos 2{pi + Ap — J/) (A6) 

D = dv 2 V 2 e{pi)e-[pi + A(())(cos2((()i - J/)) (A7) 

X (cos 2{(pi J- Ap — Jt) 

It is easy to show that the integrals of B and C are 
equal to 0 and the integral of D leads to: 


1 

TT 




2r!2U2 COS 2 A()) (A8) 

X f e{pi)e{pi + Ap) 

J —TT 


Therefore, 

Pfg{AP) = (y dpie{pi)e{pi + Ap)'j (A9) 
X (1 J- 2v2V2 cos 2Ap) 


In this section, we analytically derive the weighting fac¬ 
tor introduced in Eq. (10). We start from the azimuthal 


In a similar way one can calculate the Ap distribu¬ 
tion of mixed BG pairs produced without reaction plane 
binning: 



34 




(AlO) 
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t/2 
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xd(l)id(j)2P{(t>i — 4'i)P((/)2 — 'I'2) 

where 4 >i[ 2 ) and 4'i(2) represents the azimuthal angle of 
particle 1(2) and the reaction plane azimuthal angle of 
the events from which the particles are taken. Replacing 
P(</.-4') by (ig: 


Pm/x(A 0 ) 

= ^ / ddfi d4'2 / 

^ 7r/2 77/2 

xd(l)idME + P + G + H) 


E = e(^i)e{(/>i + A(j)) 


(All) 


(A12) 


F = 2u2e((/>i)e((^i + A(p) cos2((^i — 4'i) (A13) 


G = 2w2e(?f>i)e((/'i + A(j)) cos 2{(j)i + Acj) — 4'2)(A14) 


where {CB^ _/++/_) are the unlike-sign and like-sign 

integral yields or average numbers of pairs per event. 

We start from the case without elliptic flow. Then, 
as proven in Ref [23], if e’*' and e~ are always produced 
in pairs independent of each other, the average number 
of unlike-sign and like-sign combinatorial pairs can be 
calculated as: 

(CR-i—) = [sp + e+(l — £p)][£p -|- e_(l — £p)] (B2) 
x({N^)-{N)) 

{CB++) = + e+(l - ep)]\{N^) - (N)) (B3) 

(CP__) = + e-(l - ep)]"((A") - (N)) (B4) 

where Sp is the probability to reconstruct both tracks of a 
pair, £+/- is the probability to reconstruct only a single 
track and N is the number of pairs in an event. 


If Sp/j^j- are assumed to b e const ants, Eq. (Bl) can 
easily be proven from Eqs. ( B2|B4 ). However, in the 
presence of elliptic flow, the probabilities depend 

on the reaction plane angle: 

£p/+/_('0)=y d(()£p/+/_((/))(l + 2u2COs(((i-'(/')) (B5) 


H = 4u2V2e((^i)e((()i + Acj)) cos 2((^i — d^i) (A15) 
X cos + A(f) — 4>2) 

Because F, G and H are again easily proved to be 0, 
Pmix{A4>) can now be written as: 

PMix{A(j)) = f d4>ie{(l)i)e{(j)i + A(l)) (A16) 

J — TT 

The weighting factor to introduce the flow correlation 
into the mixed BG pairs is then given by the ratio be¬ 
tween Eq. (AfO) and Eq. (A16): 


w{A(j)) = 


PpciAcj)) 


PMixiA4>) 

= f -I- 2 v 2 V 2 cos 2Acj) 


(A17) 


Appendix B: Violation of CB^ _= 2y/CB++CB-- due 

to flow 


In this appendix, we demonstrate that the combina¬ 
tion of elliptic flow and nonuniform detection efficiency 
violates the well-known relation between unlike-sign and 
like-sign combinatorial background: 


{CB+_) = 2^{GB++){CB__) 


{GB+_i^)) = [Ai)l))B{i,)] X i{N^) - (TV)) (B6) 

{CB++m = ^[Ai)l))r X {{N^) - (TV)) (B7) 

{GB..{cP)) = liBicP)]^ X ((TV2) - (TV)) (B8) 

Ajc/j) = EpO)) + e+i-ijjjil - EpO))) (B9) 

R(i/>) =£p (!/')+£-(■!/') (l-£p(V')) (BIO) 


Taking the average over tjj within [—f, f] gives: 

{GB+_) = ((TV^) - (TV)) J di) A{i))B{i)) (BlI) 
(CR++) = ^(TV^) - (TV)) I d^ A{)l)f (B12) 

(CR__) = ^(TV^) - (TV)) I di) B{)j) f (B13) 

Using the Cauchy-Schwarz inequality, one obtains: 

< J d^P Ajijjf (B14) 


dij) Ajip)B{)/)) 


dtf) Bjif))"^ 


(Bl) 


and consequently, 

{GB+_) < 2,/{GB++){GB__) 


(B15) 
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Appendix C: A second, independent analysis 


A subset of the data, 4.8 x 10® MB events, was ana¬ 
lyzed by a second independent team. The second analysis 
follows the analysis strategy presented in Ref. [23], but 
includes the information provided by the HBD and other 
important improvements developed in this work. 

In this appendix we present the key features of the 
second analysis with an emphasis on the most important 
differences to the main analysis: (i) the HBD underlying 
event subtraction and cluster algorithm, (ii) the electron 
identification cuts and (iii) the background normaliza¬ 
tion. All analysis steps not explicitly mentioned are iden¬ 
tical between the two analyses. In particular, identical 
cuts on the acceptance and inactive detector areas, and 
the same pair cuts are applied. At the end of this ap¬ 
pendix we discuss the efficiency correction and compare 
the results of both analyses. 

The net number of photo electrons in an HBD cluster 
was calculated with a different algorithm than discussed 
in Section |H D[ using a local estimate of the scintillation 
background rather than a module average. As an electron 
typically fires three HBD readout cells, 3-cell triplets are 
used to initiate the cluster search. All possible triplets 
are formed. The photo-electron background due to scin¬ 
tillation light is estimated by the median amplitude in 
the first and second neighboring cells around the triplet. 
The background subtracted triplet charge is calculated 
as: 

qnet = qt-Ax^^£^^^ (Cl) 

where qt is the total charge in the triplet. At the number 
of cells with charge in the triplet, and (<?/«), (<?sra) are the 
median charge in the first and second neighboring cells, 
respectively. Only triplets with 0 < q^t < 60 p.e. are 
recorded. 

Electron candidates are projected to the HBD, and 
triplets within 1.5 cm of the track are merged to form 
a cluster. The net charge of the cluster qr- is calculated 
starting from the sum of the charge of all cells in the 
cluster: 

n — n —A V (^/") + (qsn) 

Qr — Qtotclust -^clust ^ 2 


where qtotclust is the sum of the charge of all cells in the 
cluster, Adust is the number of cells in the cluster, {qfn), 
(qsn) are again the median charge per cell in the first and 
second neighbors but now around the cluster. 

This analysis uses a number of sequential one¬ 
dimensional cuts to identify electrons. The variables 
used for the electron identification are dehned in Sec¬ 
tion jlTlC^ The following cuts are used: 


• nO > 2: The exclusion of RICH photo-multipliers 
fired by background electrons (Section HIC2) is 
not used in this analysis. 


• chi2/npe0 < 20 

• emcsdr < 3 

• |dep| < 2 

• ’^TOF <l-5cr: Calculated based on the time-of- 
flight measured by either the EMCal or the TOE-E 
detectors. 

• 10 < ( 7 r < 40 p.e.: Cluster charge as defined in Eq. 

With these cuts, a purity of the electron sample of 
86 % is achieved for the most central bin, which quickly 
increases to above 99% for the most peripheral collisions. 

The combinatorial background is calculated by event 
mixing. We use the method outlined in |23|j but in¬ 
cluded the weighting for the azimuthal anisotropy as im¬ 
plemented in the main analysis and described in Section 
m E 2| Eor the correlated background both analyses use 
the same MC simulations. Eor cross-pairs and jet-pairs 
the simulated pairs were reanalyzed with the track selec¬ 
tion cuts and HBD cluster algorithm mentioned above. 
The shapes of the mass spectra are consistent within sys¬ 
tematic uncertainties for the two analysis methods. For 
the electron-hadron and BB contributions the simulated 
pairs were not reanalyzed. 



FIG. 37. (Color online) The ratio of the foreground like-sign 
pairs to the sum of combinatorial and correlated pair sources 
in centrality bins 0%-10%, 10%-20%, 20%-40%, 40%-60% 
and 60%-92%. 

The normalizations of all the background components 
were fitted simultaneously to the full mass and px range 
of the like-sign spectra: 

-= qqBG^^ - -l- oiCHf-i- (C3) 


• disp < 5.5 cm 
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FIG. 38. (Color online) Comparison of final spectra from the 
main (M) and second (S) analyses. 


+02 + 03 

- + - 

The parameters ai are the individual normalization con¬ 
stants. Figure [37| shows the like-sign foreground divided 
by the sum of all background sources for the five central¬ 
ity classes. The uncertainty on the combinatorial back¬ 
ground normalization is shown as a gray band on each 
panel. No systematic deviation from unity is observed, 
indicating that the sum of the different background com¬ 
ponents gives a sufficiently accurate description over the 
mass range up to 2 GeV/with no indication of any 
shape variation within the shown uncertainties. Above 
2 GeV/c^ the statistical significance makes a comparison 
at the shown scale meaningless. 

After fixing the normalization of all background 
sources so that a satisfactory description of the like-sign 
pairs is achieved, the analysis is extended to unlike-sign 
pairs. The normalizations for the unlike-sign cross-pairs, 
jet-pairs and electron-hadron pairs are taken from Eq. 
(G4). For the combinatorial unlike-sign pairs we use 


unlike-sign mixed event pairs. The normalization is also 
taken from Eq. (G4), but needs to be corrected to ac¬ 


count for the different effect of the pair cuts on like- and 
unlike-sign pairs as done in Ref. [53]. 


To estimate the uncertainty on the raw yield due to the 
background subtraction one needs to consider the signal- 
to-background ratio S/B. The uncertainties on the 
are multiplied hy B/S and added in quadrature. This 
results in ^ 55% systematic uncertainties at 0.6 GeV/c^ 
for MB collisions. 

We factorize the efficiency into 3 terms, which are de¬ 
termined separately. 

total _ ^ TOF ^ embed (C^A\ 

^pair ^pair ’ ^pair ’ ^pair 

The first factor describes the effect of all reconstruction 
algorithms and cuts except for the time-of-flight cut and 
the centrality dependence of the reconstruction efficiency 
in the central arms, which are treated separately. It is 
obtained by a MG simulation of e+e” pairs that are pro¬ 
cessed through the full PHENIX detector simulation, in¬ 
cluding the HBD. The simulated HBD hits are embedded 
into real HBD data as discussed in Section fill FI These 
events are then analyzed with the same electron identi¬ 
fication, fiducial, and pair cuts used in the independent 
analysis, with exception of the time-of-flight cut. The 
systematic uncertainty of Cpair is about 12%. It was de¬ 
termined from the measured yield of pairs in the Dalitz 
decay region when varying electron identification cuts in 
a way that changes the raw pair yields by factors between 
0.5 and 1.5. 

The efficiency is determined from tracks mea¬ 

sured in peripheral collisions, where the hadron contam¬ 
ination is negligible, by comparing data obtained with 
a 1.5 cr cut to the case with no time-of-flight cut. We 
find that on average the TOF efficiency for tracks is 93% 
above 0.4 GeV/c, but drops to 80% at 0.2 GeV/c inde¬ 
pendent of centrality. This drop results from a failure of 
the electronics to properly record time for low amplitude 
signals. In the main analysis this issue was avoided by 
treating tracks with no time information separately. The 
systematic uncertainty due to this cut is a few percent at 
0.6 GeV/c^. 

The efficiency was determined by embedding 

MC-simulation tracks into the data of all used central 
arm detectors and analyzing these embedded tracks using 
the same cuts as used in the data. The values are found to 
be very similar to those derived in the main analysis. For 
central collisions an additional 8% systematic uncertainty 
is added. 

Gompared to the main analysis, the total reconstruc¬ 
tion efficiency is a factor of ~2 smaller for central 
collisions. The difference drops to ~30% for the most 
peripheral collisions. 

The fully corrected mass spectra from the independent 
analysis are compared to those from the main analysis 
in Fig. for all five centrality bins. The results are 
consistent within uncertainties. 
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